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Abstract 



The adiabatic approximation is a natural approach for the description of phenomena 
induced by low frequency laser radiation because the ratio of the laser frequency to the 
characteristic frequency of an atom or a molecule is a small parameter. Since the main 
aim of this work is the study of ionization phenomena, the version of the adiabatic ap- 
proximation that can account for the transition from a bound state to the continuum must 
be employed. Despite much work in this topic, a universally accepted adiabatic approach 
of bound-free transitions is lacking. Hence, based on Savichev's modified adiabatic ap- 
proximation [Sov. Phys. JETP 73, 803 (1991)], we first of all derive the most convenient 
form of the adiabatic approximation for the problems at hand. Connections of the ob- 
tained result with the quasiclassical approximation and other previous investigations are 
discussed. Then, such an adiabatic approximation is applied to single-electron ionization 
and non-sequential double ionization of atoms in a strong low frequency laser field. 

The momentum distribution of photoelectrons induced by single-electron ionization is 
obtained analytically without any assumptions on the momentum of the electrons. Previous 
known results are derived as special cases of this general momentum distribution. 

The correlated momentum distribution of two-electrons due to non-sequential double 
ionization of atoms is calculated semi-analytically. We focus on the deeply quantum regime 
- the below intensity threshold regime, where the energy of the active electron driven by 
the laser field is insufficient to coUisionally ionize the parent ion, and the assistance of 
the laser field is required to create a doubly charged ion. A special attention is paid 
to the role of Coulomb interactions in the process. The signatures of electron-electron 
repulsion, electron-core attraction, and electron-laser interaction are identified. The results 
are compared with available experimental data. 

Two-electron correlated spectra of non-sequential double ionization below intensity 
threshold are known to exhibit back-to-back scattering of the electrons, viz., the anticorre- 
lation of the electrons. Currently, the widely accepted interpretation of the anticorrelation 
is recollision-induced excitation of the ion plus subsequent field ionization of the second 
electron. We argue that there exists another mechanism, namely simultaneous electron 
emission, when the time of return of the rescattered electron is equal to the time of liber- 
ation of the bounded electron (the ion has no time for excitation), that can also explain 
the anticorrelation of the electrons in the deep below intensity threshold regime. 

Finally, we study single-electron molecular ionization. Based on the geometrical ap- 
proach to tunnelhng by P. D. Hislop and I. M. Sigal [Memoir. AMS 78, No. 399 (1989)], 
we introduce the concept of a leading tunnelling trajectory. It is then proven that leading 
tunnelling trajectories for single active electron models of molecular tunnelling ionization 
(i.e., theories where a molecular potential is modelled by a single-electron multi-centre 
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potential) are linear in the case of short range interactions and "almost" linear in the case 
of long range interactions. The results are presented on both the formal and physically 
intuitive levels. Physical implications of the proven statements are discussed. 
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Chapter 1 
Introduction 



In this thesis we explore one- and two-electron ionization phenomena in a strong laser 
field. Since the ratio iojujat ~ 0.1, where uj is, say, the frequency of Ti:sapphire lasers - 
most commonly used tabletop laser systems in laboratories around the world jl], and lOat 
is the characteristic atomic frequency, the adiabatic approximation ought to be a perfectly 
suitable theoretical tool for the description of phenomena induced by such low-frequency 
laser fields. 

Roughly speaking, the adiabatic approximation can be introduced as follows. Once 
the frequency of the external laser field is much lower than the characteristic atomic fre- 
quency, uj <^ ujat-, an approximate solution of the Schrodinger equation can be found by 
means of averaging over atomic (internal) degrees of freedom. Therefore, the adiabatic 
approximation is a method of constructing an asymptotic expansion of the solution of the 
non-stationary Schrodinger equation in terms of the small parameter u/uat- 

Born and Fock [2] founded the theory of the adiabatic approximation for a discrete 
spectrum by formulating the adiabatic theorem. Landau [3l H] estimated the probability of 
nonadiabatic transitions between discreet states. However, the leading-order asymptotic 
result for such a quantity was obtained by Dykhne f5j . Afterwards, nonadiabatic transitions 
between discreet states were thoroughly analyzed by many authors [SI El E] (for reviews 
see References [9l [ini HI])- Summarizing research in this area, one may conclude that 
nonadiabatic transitions in discreet spectra are quite well studied. 

However, transitions from a discreet state to a continuum have been the subject of 
on-going investigations for many decades. Despite much work on this topic, a universally 
accepted approach is lacking. The direct generalization of the Dykhne method was devel- 
oped by Chaplik [T21 US]- Exactly solvable models were reported by Demkov and Oserov 
[Tit US], Ostrovskii [16], Ostrovskii et al. [HI HE], and Nikitin ^ [20] (for review see, 
e.g.. Reference [21]). The advanced adiabatic approach was introduced and widely em- 
ployed by Solov'ev [221 ISSl [21]. Finally, Tolstikhin recently developed a promising version 
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of the adiabatic approximation for the transitions to the continuum [25] by employing the 
Siegert-state expansion for nonstationary quantum systems [261 [271 [28] . Yet, the approach 
presented in Reference [23] is hmited to finite range potentials. 

The adiabatic approximation deserves special attention because it is one of a few known 
non-perturbative approaches in quantum mechanics. Strong field physics demands for such 
methods to describe nonlinear phenomena induced by a strong laser field. 

In Chapter [2| we derive the amplitude of nonadiabatic transitions from a bound state 
to the continuum within the Savichev modified adiabatic approximation [29]. Then, the 
properties of this amplitude are scrutinized, and a connection between the adiabatic and 
quasiclassical approximations is established. The obtained amplitude is the main theoreti- 
cal formalism of this work, which is applied to strong field one- and two-electron processes 
in subsequent chapters. 

The most common process that occurs when laser radiation is exerted on an atom is 
single-electron ionization. Although initial theoretical understanding of strong field ion- 
ization was put forth by Keldysh [30j as early as 1964, many questions remain unresolved. 
As far as single-electron ionization in the presence of a linearly polarized laser field is con- 
cerned, there are two important topics. The first, namely the ionization rate as a function 
of instantaneous laser phase, was studied in depth by Yudin and Ivanov [3T] (see also Ref- 
erences [32l [33]). However, their result assumes zero initial momentum of the liberated 
electron. Effects due to nonzero initial momentum have yet to be included. The second 
topic pertains to the single-electron spectra, that is, the ionization rate as a function of the 
final momentum of the electron. Despite much work on this topic (see, for example. Refer- 
ences [311 [351 [3S] and references therein) a universal formula is absent and the discussion 
is still ongoing. Among the most accurate results, Goreslavskii et al. [36] have obtained 
an expression for the complete single-electron ionization spectrum, but without consider- 
ation of the laser phase. In Chapter |3| we derive a more general formula that includes the 
dependences on both the instantaneous laser phase and the final electron momentum. 

In strong low-frequency laser fields, following one-electron ionization of an atom or a 
molecule, the liberated electron can recoUide with the parent ion [371 [38] . The electron acts 
as an "atomic antenna" [37], absorbing the energy from the laser field between ionization 
and recoUision and depositing it into the parent ion. Inelastic scattering on the parent ion 
results in further collisional excitation and/or ionization. Liberation of the second electron 
during the recollision - the laser-induced e-2e process - is known as nonsequential double 
ionization (NSDI). 

The phenomenon of NSDI was experimentally discovered by Suran and Zapesochny 
[39] for alkaline-earth atoms (for further experimental investigations of NSDI for alkaline- 
earth atoms, see, e.g.. References [101 [HI [121 [131 [lH [IS] ) • In this case, autoionizing double 
excitations below the second ionization threshold were shown to be extremely important. 
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For a theoretical study of these effects, see, e.g.. Reference |16]. For noble gas atoms, 
nonsequential double ionization was first observed by L'Huillier et al. (see, e.g.. References 
[l?! HH]). The interest to the phenomenon of NSDI grew rapidly after it was rediscovered in 
1993-1994 [iniEn]- Recently, correlated multiple ionization has also been observed [^152]. 
The renewed interest in NSDI has been enhanced by the availability of new experimental 
techniques that allow one to perform accurate measurements of the angle- and energy- 
resolved spectra of the photoelectrons, in coincidence. Such measurements play a crucial 
role in elucidating the physical mechanisms of the NSDI. 

From the theoretical perspective, direct ah initio simulations of the photoelectron spec- 
tra corresponding to NSDI in intense low-frequency laser fields represent a major challenge. 
Only now such benchmark simulations have become possible [53] for the typical experimen- 
tal conditions (the helium atom, laser intensity I ~ 10^^ W/cm^, laser wavelength A = 800 
nm). In addition to these calculations, tremendous insight into the physics of the problem 
has been obtained from classical simulations performed in References [Sl f I55ll56ll57t l58 | 159] . 
These papers have demonstrated a variety of the regimes of nonsequential double and triple 
ionization. Not only do these simulations reproduce key features observed in the experi- 
ment, they also give a clear view of the (classical) interplay between the two electrons, the 
potentials of the laser field, and of the ionic core. They also show how different types of 
the correlated motion of the two electrons contribute to different parts of the correlated 
two-electron spectra. 

The physics of double ionization is different for different intensity regimes, separated 
by the ratio of the energy of the recoUiding electron to the binding (or excitation) energy 
of the second electron, bound in the ion. 

According to classical considerations, the maximum energy which the recolliding elec- 
tron can acquire from the laser field is ~ 3.2Up [38j, where Up = {F/2ujY, F is the 
laser field strength, and uj is the laser frequency (unless stated otherwise atomic units, 
h = nif. = \e\ = 1, are used throughout the work). Hence, NSDI can be divided into two 
types: if the intensity of the driving laser field is such that the recolliding electron gains 
enough kinetic energy to coUisionally ionize the parent ion - the above intensity threshold 
regime, and when such kinetic energy is insufficient to directly ionize the ion - the below 
intensity threshold (BIT) regime. The former regime is thoroughly studied experimentally 
as well as theoretically (see, e.g.. References jnOlEIlEaESlElIMlESlEZllSilMlEniin^ 
and references therein). 

NSDI BIT, being a most challenging regime, is currently of an active experimental 
interest [721 EH ESI EH E21 ESI ES] • In this regime, existing classical and quantum analysis 
(see, e.g.. References [SSlEZlEZlESlEniEniEIllSSEB]) demonstrates two possibilities 
of electron ejection after the recoUision. First, the two electrons can be ejected with 
little time delay compared to the quarter-cycle of the driving field. Second, the time delay 
between the ejection of the first and the second electron can approach or exceed the quarter- 
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cycle of the driving field. In these two cases, the electrons appear in different quadrants 
of the correlated spectrum. If, following the recoUision, the electrons are ejected nearly 
simultaneously, their parallel momenta have equal signs, and both electrons are driven by 
the laser field in the same direction toward the detector. If, following the recoUision, the 
electrons are ejected with a substantial delay (quarter-cycle or more), they end up going 
in the opposite directions exhibiting the phenomenon of anticorrelation of the electrons. 
Thus, these two types of dynamics leave distinctly different traces in the correlated spectra. 

In Chapters [4], we develop a fully quantum, analytical treatment of NSDI BIT. It is 
important that our approach takes into account all relevant interactions - those with the 
laser field, the ion, and between the electrons - nonperturbatively. The case in which 
the two electrons are ejected simultaneously, i.e., the process of simultaneous electron 
emission (SEE), is considered. We show that in this case the correlated spectra bear clear 
signatures of the electron-electron and electron-ion interactions after ionization, including 
the interplay of these interactions. These signatures are identified. In agreement with 
previous studies, the mechanism of SEE manifests itself in the correlation of the electrons 
- the electrons are moving in the same direction after NSDI. 

However, in Chapters [5} we demonstrate that if the intensity of the laser field is lowered 
such that we enter the deep BIT regime of NSDI, SEE can be responsible for the anticorre- 
lation of the electrons. This novel mechanism is alternative to the widely accepted point of 
view that the anticorrelation of the electrons are caused by recollision-induced excitation 
of the ion plus subsequent field ionization of the second electron (RESI). Nevertheless, 
SEE and RESI are by no means mutually exclusive processes; they both contribute to the 
complex and diverse phenomenon of NSDI BIT. 

Recent advances in experimental investigations of single-electron molecular ionization 
in a low frequency strong laser field [HI [851 ESI ETJ EHJ [89l [901 IHI] have created a demand 
for a theory of this phenomenon [n21[n31[nil[nSl[nni[SZl[nHl[nSl [M [Ml [102] • As far 
as low frequency laser radiation is concerned, one can ignore the time-dependence of the 
laser and consider the corresponding static picture, which is obtained as a; — )■ 0. In 
this limit, single electron molecular ionization is realized by quantum tunnelling. This 
approximation is valid from qualitative and quantitative points of view, and it tremendously 
simplifies the theoretical analysis of the problem at hand. Such single active electron 
approaches to molecular ionization, where an electron is assumed to interact with multiple 
centres that model the molecule and a static field that models the laser, are among most 
popular. Analytical and semi-analytical versions of these methods, which are based on the 
quasiclassical approximation [931 ESI [Ml [991 HOP] are indeed quite successful in interpreting 
and explaining available experimental data. However, these quasiclassical theories heavily 
relay on the assumption that the electron tunnels along a straight trajectory. Despite its 
wide use, the reliability of this conjecture has not been verified. 

In Chapter |6| we study the reliability of this hypothesis. Relying on the geometrical 
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approach to many-dimensional tunnelling by Hislop and Sigal |103[ 11041 11051 1106] , which 
is a mathematically rigorous reformulation of the instanton method, we first introduce the 
notion of leading tunnelling trajectories. Then, we analyze their shapes in the context of 
single active electron molecular tunnelling. It will be rigorously proven that the assumption 
of "almost" linearity of leading tunnelling trajectories is satisfied in almost all the situations 
of practical interest. Such results justify the above mentioned models and open new ways 
of further development of quasiclassical approaches to molecular ionization. 
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Chapter 2 

Adiabatic Approximation 



2.1 General Discussion 



Our investigations are based on a seminal result that ought to be summarized foremost. 
Following the Solov'ev advanced adiabatic approach [22l [231 El] , Savichev p9] proved the 
following. If the adiabatic state \ipi{ip)) and the corresponding adiabatic term Ei{ip), 



are known, then the solution {"^{t)) of the nonstationary Schrodinger equation 

td\^{t)) /dt = H{^)\^{t)), 
where ip = cut is a. phase of the laser field, subjected to the initial condition 



I— > — oo 

has the following form within the adiabatic approximation {u <^ 1) 

1 



2nuj 
exp 



dEdif' \^P^i^^{E))) X 



(jJ 



Eip' 



Ei{(fi)d(fi - iEt 



[1 + OM]. 



(2.1) 



(2.2) 



(2.3) 



(2.4) 



Here ^Pi{E) is the analytical continuation of the inverse function of Ei{(p). Not e th at no 
assumption^ on a form of the Hamiltonian H are required to derive Equation (2.4) from 
the nonstationary Schrodinger equation (2.2). 



^ Besides some tacit requirements such as the analyticity of both the adiabatic terms and the states. 
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The integral over E in Equation (2.4) can be interpreted as a generalization of the 



Born-Fock expansion to the case of a system with a continuos spectrum. The integral 



over (f' in Equation (2.4) can be calculated by the saddle point approximation without 



changing the accuracy of Equation (2.4); however, the original integral representation is 



more advantageous and should be left unaltered. 

For the sake of completeness and clarity, we shall present the derivation of Equation 



(2.4). Let us seek the solution of the Schrodinger equation (2.2) in the form 



|*(t)) = j dEg{E)exp{-tEt) \U^^{E))) 
where the unknown function g{E) can be represented as 



(2.5) 



giE) = ^ I dt expitEt)Git,E), G{t,E) = {UVr{E)) \^{t)) 



(2.6) 



where the normalization condition \ijji{ip)) = 1 is assumed for almost all ip. To 

derive the equation for G{t,E), we recall that |^(t)) satisfies the Schrodinger equation 



(2.2); hence. 



{Ucp,iE))\^^^\^it)) = {U^,iE))\Hi^)\^it)) 
^^G{t,E) = {^Mim^§-^\m)■ 



(2.7) 



Having substituted Equations (2.5) and (2.6) into the system of equations (2.7), we get 



j j dE'dt' eM^E'{t' - t)]G(t', E')E'{,p,{^,{E)) \^,{^,{E'))) 
= jj dE'dt' eM^E'{t' -t)]G{t',E') (^,((^,(E))|i^((^)|^,((^i(E'))), 
^^G(t,E) = ^jj dE'dt' eMiE'it' - t)]G{t' , E')E' {iJi{^,{E)) \^P^{^^{E'))) 
Introducing the variables ip = ut and ip' = out', we have 



JJ dE'dif' exp 
= JJ dE'd^' exp 
toj^^Gi^,E) 



-E'{^' - ^) 



G{v\E')E'{^lj,{v,{E))\Uv,{E'))) 



-E'{^' - ^) 



G{v\E') {Uv^m\ H{V) mv^{E'))) , 



2ti 



dE'dif' exp 



-E'{ip' - ^) 



G{^',E')E'{i;,{^,{E))\i^M{E'))) 



Substituting the following ansatz into the equations above 



G{^,E) = [Xo {E, ^) + LoXi (E, yp) + ■ ■ ■ ] exp 



I 



Q{ip')dip' 



(2.8) 



then performing integration over E' and t' by means of the saddle point approximation 
and collecting terms in front of the zeroth power of u, we obtain 

Xo{E,^) = Xo{Q{v>),V>)mV>^iE)) \Uv^iQiv))))■ (2.9) 

Whence, 

Q{V)=E,{^), Xo{E,<f) = \ i'i{^^{E,{^)))) = {i,,{^,{E))\i,,{^)) . (2.10) 



Substituting Equations (|2.8|) and (|2.10|) into Equation (|2.6|), we have 
1 



5(^) 



27ra; 



E^ 



Ei{f')d^' 



[l + Oiu)]. (2.11) 



According to the saddle point approximation, the only neighbourhoods of importance in 
the integral (2.11 ) are those where the derivative of the exponent with respect to vanishes. 
In these neighbourhoods the matrix element {ilJi{ipi{E)) \^|Ji{ip)) = 1; therefore, the integral 
(2.11) is equivalent, up to a term of order of u, to the expression 

1 



9{E) 



271UJ 



dip exp 



E^ 



E,{^')d^' 



(2.12) 



Recalling Equation (2.5), we conclude that Equation (2.4) is finally derived. Note that a 
proper choice of complex integration paths for the integrals in Equation ( 2.4[ ) will ensure 
that the wave function (2.4) indeed satisfies the initial condition (2.3). (A more detailed 
version of the derivation of Equation (2.4) presented above can be found in Reference j29j . ) 



However, one must be cautious regarding Equation (2.4). As far as rigorous asymptotic 
analysis is concerned, it is incorrect to assume that the wave function (2.4) obeys the 
Schrodinger equation (2.2) even though Equation (2.4) was obtained from Equation (2.2). 
On the contrary, the validity of the solution (2.4) must be verified independently. 

Employing the equality 



H{^i{E)) m^.{E))) = E,{^,{E)) m^i{E))) = E \U^,{E))) 
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and substituting Equation (2.4) into Equation (2.2), one readily shows that 
.d 



\^it)) 



1 



27ra; 
exp 
1 



2 



ff 



27TU J J 

exp 



Elf' - j Ei{if)difj - iEt 
dEdip' H{ip) X 



[1 + OH], (2.13) 



Elp 



Ei{Lp)dLp - iEt 



[l + 0{u)]. (2.14) 



Hence, in order to prove that the wave function (2.4) indeed satisfies Equation (2.2), we 
ought to demonstrate that H{(p) \ipi{(p')) ~ H{ip') \i/ji{ip')) . 

Using the analyticity of H{f) and {ipiif)), we write 



n=0 



Since 



d'^Hiip) 
dip'^ 



d 


( d^'-^Hip) 






d 




f d''-^H{p)) 


dip 







d d^'-^Hiip) 
dip'' dip^'^ 



d'^-^Hiip) f d 



dip"-~^ \dip 



we obtain 

H{ip + e)\^lj,iip)) 



OO 

y- 

n=0 



d 

dip ' 



d 

dip ' 



V 

n times 



d 



dp 



exp e— H{ip)exp -e— \ipi{ip)) . 



d 



dip 



(2.15) 



Whence, we lamentably observe that the wave function (2.4) does not obey the Schrodinger 

theless, if 

Hip)] = 0(1), (2.16) 



equation (2.2) in a general case. Nevertheless, if 

d 



dip ' 
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then 



H{^ + e)\U^))=H{^) |^,(^))+o(l) 



(2.17) 



i.e., Equation (2.16) is the assumption on the form of the Hamihonian that guarantees 
that the wave function 



exp 



Ei{Lp)dip - lEt 



[1 + 0(1)] 



(2.18) 



is indeed a solution of the Schrodinger equation (2.2). [Note a minor difference between 



Equations (2.4) and (2.18). 



What are the imphcations of the condition (2.16) for strong field physics? Let us specify 



the Hamiltonian. A typical hamiltonian of a system interacting with an external laser field 
reads (in the length gauge) H{^) = p^/2 + V{r) + r • F(c<;, ip), where F(co', ip) denotes the 
laser pulse. Then, the following 



lim F(ci;, ip) 



lim dF(uj, 'p)/dip = 0, 



(2.19) 



would satisfy Equation (2.16). A strictly periodic laser pulse, such as F{u,(p) = F^sinyj, 
does not satisfy the condition (2.19). However, a pulse with a Gau ssian envelope, e.g., 
F{uj,ip) = F^exp [—{^)'^] sirup, obeys it. Therefore, the condition (2.16) demands that 



the laser field used has to have an envelope, which is a realistic and, perhaps, even tacit 
requirement. 

Due to the connection (see References [22]) between the Savichev adiabatic approach 



and the Solov'ev advanced adiabatic method, the condition (2.19) also applies to the latter 
method. 

Finally, it is important to note that results obtained within different versions of the 
adiabatic approximation are in fact the same. This follows from the uniqueness of the 
asymptotic expansion (see, e.g.. Reference |107] ). Furthermore, the results are also gauge 
independent. Hence, the choice of the gauge as well as the choice of the specific adiabatic 
method is merely the issue of convenience. 



2.2 The Derivation of the Amphtude of Non-adiabatic 
Transitions 



Let \i) and |/) be stationary states (for specification see Equation (2.24) and the comment 
after), and we shall assume that the quantum system with the Hamiltonian H is in the 
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state l^) at t = — oo. The main aim of this section is to obtain the general form of the 
transition amphtude OJtj^./ that the given quantum system will be found in the state |/) 
at t = +00. 

Before going further, we are to introduce notations. First, we arbitrarily partition the 
Hamiltonian H: 

H{ip)=H,{ip) + V{ip). (2.20) 



Second, we denote by |\E'ij(t)) the solutions of Equation (2.2) such that |\E'j(— 00)) = 



l%(+oo)) 
equation 



I/); similarly, |<l>jj(t)) are the solutions of the nonstationary Schrodinger 

with the initial conditions: |$j(— 00)) = \i) and |<l'/(+oo)) = |/), correspondingly. 

Having defined all necessary functions, we introduce two equivalent forms of the tran- 
sition amplitude OJtj^/ by employing the corresponding version of the S*- matrix (see, e.g.. 
References [66 | I108[ 1109] ): the reversed time form (sometimes called the "prior" form) 



(r) 



[^>f{t)\V{ujt)\^,{t))dt 



(2.21) 



and the direct time form (the "post" form) 



— z 



{^j{t)\V{ujt)\^>m)dt. 



(2.22) 



It is noteworthy to recall the physical interpretation of Equations (2.21) and (2.22). The 
terms (^1^/(^)1 V{ujt) |$i(t)) and ($/(t)| V{ujt) \^i{t)) can be regarded as the amplitudes of 
quantum "jumps," which occur at the time moment t. The integrals over t convey that 
these jumps take place at any time. 

Introducing the adiabatic state |0j (</))) and term Ef{(f) of the Hamiltonian Hq, 

Ho{cp)\M^)) = Eficp)\Mv)), (2.23) 



the wave function |$j(t)) can be readily presented in the form of Equation (2.18). In further 
investigations, we employ the post form [Equation (2.22)], and thus we shall assume that 



|^,(-oo))^|z), \4>f{+oo))^\f) 



(2.24) 



In the case of the prior form [Equation (2.21)], condition (2.24) has to be substituted by 
|0j(— 00)) = \i) and |'?/'/(+oo)) = |/), where \(j)i{(p)) and li^fi^)) are adiabatic states of 
the Hamiltonians Hq and H, correspondingly. 
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Substituting the asymptotic representations [Equation (2.18)] of the wave functions 



|$/(t)) and \'^i{t)) into Equation (2.22), we obtain 



{d) 



(27r)2a;3 



/(z)e*^(")/"d^z [1 + o(l)] 



(2.25) 



where z = {£,ri,£',ri',(f) is a five-dimensional vector, d^z = dSdrjdS'dri'dLp, /(z) = 

i^ji£'))\ Vi^) i^,{£))), and 5(z) = £ (r^ - ^)+£' (^ - r^')"/' ^^,.(0^^+/"' i?/(0^^e 
Bearing in mind that 1/u is a large parameter, the five-dimensional integral in Equation 



(2.25) can be calculated by means of the saddle-point approximation. Finally, the post 



form of the transition amplitude within the adiabatic approximation reads 



id) 



2tt 




(0/(<y9o)| V{^q) \lpi{(pQ)) 



dip 



X exp 



[Ef{^)-E,{^)]d^\[l + o{l)] 



(2.26) 



where denotes the summation over simple saddle points y^o, i-e., solutions of the 

equation 



d 

dip 



d 



(2.27) 
(2.28) 



The physical interpretation of the sum over ip^ is as follows: quantum jumps occur only 
at isolated time moments to = Vo/^i when the jumps are most probable; hence, to are 
called "transition times." Note that the given interpretation deviates from the physical 



meaning of the time integral in Equation (2.22). 



Some general remarks on Equation (2.26) are to be made: 



i. 9?o is usually a complex solution of Equation (2.27); therefore, saddle points (po with 



negative imaginary parts should be ignored because such points make exponentially 
large contributions to the amplitude, which leads to unphysical probabilities. 



ii. If (y9Q, (y9Q, . . . , v^q a^c solutious of Equation (2.27), such that Im {ip>\) > Im {ipD > 



. . . > Im [iPq] > 0, then all but the single term that corresponds to the saddle point 



may be neglected in the sum over ipo in Equation (2.26). One can do so since this 
saddle point has the largest contribution to the transition amplitude. 

iii. On the one hand, the explicit form of Ef{ip) is solely determined by partitioning 
[Equation (2.20)]; on the other hand, Ei{(p) is unique for a given quantum system. 
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IV. 



V. 



The exponential factor of Equation (2.26) is similar to the exponential factor in the 
Dykhne approach O [121 [131 |6] (see also References |110l llllj ) - the methods for 
calculating the amplitude of bound-bound transitions within the adiabatic approxi- 



mation. Hence, Equation (2.26) may be considered as a generalization of the Dykhne 



formula for bound-free transitions. 

By employing an appropriate version of the saddle-point method, one can in principle 



generalize Equation (2.26) for the case when the condition (2.28) is violated. 



2.3 A Connection between the Quasiclassical and Adi- 
abatic Approximations 



Now, the connection between the amplitude [Equation (2.26)] and the method of complex 
trajectories is to be manifested. According to the method of complex quantum trajectories 
(see, e.g.. References [3l HI I110[ I112j . and the imaginary time method |113] ). to calculate 
the probability of the transition from the initial state to the final, one should first solve 
the corresponding classical equations of motion and find the "path" of such a transition. 
However, this path is complex; in particular, the transition point fq and transition time 
to at which the transition occurs are complex. Parameters Tq and to are determined by 
the classical conservation laws as shown below in this section. Next, one has to obtain the 
classical action Sf{rf,tf; Fq, to) + Si{rQ, to; fi, tj) for the motion of the system in the initial 
state from the initial position Fj at time tj to the transition point Fq at time to and then 
in the final state from fo at to to the final position fj at time t/. Finally, the probability 
of the transition is given by 



r oc exp {-2Im [S'/(F/,t/; Fo,to) + S'i(Fo,to; ri,U)]} . 



(2.29) 



Equations (2.26) and (2.29) must coincide in some region of parameters. The method of 



complex trajectories can be derived as the quasiclassical approximation of the transition 



amplitude [Equation (2.21) or Equation (2.22)]; we outline this derivation below. There- 
fore, it would be of methodological interest to establish an explicit connection between 



Equations (2.26) and (2.29). 



Without loss of generality, assuming that V{ut) is a non-differential operator, we obtain 



the quasiclassical approximation to Equation (2.22) 



(d) 



dt I dhdhfd-'Yi {f\Yf)V{ujt,Y) {y,\i) 



xFfFiexp {i [Sf{vf,tf; F,t) + Si{v,t; Fj,ti)]} 



(2.30) 
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where tf^i = ±00, Ff exp{iSf), and Fiexp{iSi) are the quasiclassical versions of the propa- 
gators with the Hamiltonian Hq and H, correspondingly. We recall that the general form 
of the quasiclassical propagator is given by 



J2 ^^"^ (r, t; r', t') exp [iS^"'^ (r, t; r', t')] , 



(2.31) 



where the sum denotes the summation over classical paths that connect the initial (r', t') 
and final (r,t) points. Therefore, usage of this form of the quasiclassical propagator, 
Fexp{iS), is justified if we assume that there is only one such path; indeed, this is the 
case in the majority of practical calculations, and thus we shall accept this assumption 
hereinafter. 



0, 
0. 



In order to reach Equation (2.29) from Equation (2.30), one has to calculate the integrals 
over r and t in Equation (2.30 ) by means of the saddle-point approximation. The equations 
for the saddle points tq and to, i.e., the transition points, read 

d [S,(r, t; v,,U) - Sf{r, t; Vf, tf)] /dt\^^,^^^^^^ 
Vr [Si{v, t; Vi, ti) - 5/(r, t; rj, tj%^^^^^^^^ 

Recalling the Hamilton- Jacobi equation 

dSijir, t; r,j, tij)/dt = - J^,/(r, j, t), 

where J^ij are classical Hamiltonians and are classical canonical momenta 

^i,/(r, t) = Vr^i j(r, t; r^j, tij), 



(2.32) 
(2.33) 

(2.34) 



we rewrite Equations (2.32) and (2.33) as the law of conservation of canonical momentum 
and the law of conservation of energy: 



^/(ro,^o) 
Jf^(ro,^/(ro,to),to) 



^(ro, ^i(ro,to),to)- 



(2.35) 
(2.36) 



Having introduced all the necessary quantities, we demonstrate the correspondence 
between Equations (2.26) and (2.29) within an exponential accuracy. Performing a simple 
transformation and using Equation (2.34), we reach 



'S'/(i"/7 tf] i"07 ^o) + Si{rQ, to; fi, ti 
•^/("^Z' ^f'^ ^i) + 'S'j(ro, tf, Fj, ti) 
•^/("^Z' ^/i "^05 ti) + Si{rQ, ti] Ti, ti) 

"to 



to 



d_ 



Si{^G, t; Tj, ti 



d_ 



Sf{ro,T;rf,tf) 



dr 



rto 

+ / [^/(ro, ^/(ro, r), r) - ^(fq, ^i(ro, r), r)] dr. 



(2.37) 
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Usually in the case of multiplioton ionization, to is complex and ro is real; moreover, 
energies along the trajectories are always real - even under the barrier. This suggests that 



the first two terms of the right hand side of Equation (2.37) affect only the phase and do 
not contribute to the probability. 

We point out that sometimes it is useful to employ a mixed representation, such as 
'^/('^Z' "^O; ^o) + Si{rQ, to; Tj, ti) 



if,(r,r,r) 



d_ 



5/(ro,r; r/,t^ 



dr 



Jti 

= Sf{rf,tf,ro,U)+ / [if,(r, r, r) + jr/(ro,^/(ro, r),r)]rfr, (2.38) 

Jti 

where is the classical Lagrangian, Si = J ^idr. 

Finally, since Jff and J^fi are the quasiclassical limits oi Ef and Ei (this will be demon- 
strated below), we conclude that the exponential factors of Equations ( |2.26 ) and (2.29) 
indeed coincide within the quasiclassical approximation. 

The wave function 



vl>,,(r,t) = j F,(r,t;r',t,)e^^'('''*^'''*'V^n(r)rfV, 



(2.39) 



is the (leading-order term) quasiclassical solution of Equation (2.2) with the initial con- 
dition \l/gc(r,ti) = 0m(i")- Employing Equation (2.34) and bearing in mind that = 
,t) does not depend on r', we obtain 



H^,, = ihdt^.c [1 + 0{h)] = .m<ic [1 + 0{K)\ . 



(2.40) 



Since we have freedom of choosing the initial condition 0j„(r), there are in general infinitely 
many wave functions [Equation (2.39)] that satisfy Equation ( |2.40 ). Comparing Equations 
(2.40) and (2.1), and taking into account the latter, we formulate the following property 



of the adiabatic term and state of a given quantum system within the quasiclassical limit: 
there exists only one adiabatic term, which is equal to the classical Hamiltonian, and any 
solution of the corresponding nonstationary Schrodinger equation is also an adiabatic state 
that corresponds to this adiabatic term (i.e., the adiabatic term is infinitely degenerate). 
Note that this property is completely ruled out once the general form of the quasiclassical 



propagator (2.31) is considered 



The property stated above, nevertheless, merely accentuates the fundamental difference 
between the quasiclassical and adiabatic approximations. As mentioned in Chapter [T| the 
adiabatic approximation allows us to obtain the solution of the nonstationary Schrodinger 
equation as an asymptotic series in terms of the small parameter u/uat'-, however, the 
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quasiclassical approximation is a method of obtaining an asymptotic expansion of the 
solution with respect to the small parameter h. These two series may be dissimilar in a 
general case. 



2.4 Summary 

Since the adiabatic approximation shall be used in subsequent chapters. It is convenient 
to conclude this section with a short summary of the main equations. 

Let the Hamiltonian of a system H{ujt) be a slowly varying function of time t, i.e., 
uj <^ 1, then the rate of a non-adiabatic transition is given by 

r^exp(^-hm p[Ef{ip)-Ei{^)]di^, (2.41) 

where ip = cut, ipo is the complex solution of the equation 

Eiiipo) = Ef{ipo) (2.42) 

with the smallest positive imaginary part, and Eij{ip) being the adiabatic terms of the 
system "before" and "after" the non-adiabatic transition. If the quasiclassical approxima- 
tion to r is sufficient, then instead of using the adiabatic terms, one can employ the total 
energies of the corresponding classical system. 
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Chapter 3 

Single-electron Ionization 



3.1 Main Results 



The Keldysh theory [30] was reformulated in terms of the Dykhne adiabatic approximation 
in Reference for the first time. Later, this approach was employed in References 

[migniM fTTSlfTTB] . 



Similarly, we shall apply the adiabatic approach [Equation (2.41)] to the problem of 
ionization of a single electron under the influence of a linearly polarized laser field with the 
frequency uj and the strength F. The initial and final classical energies for such a process 
are given by 



(3.1) 



where Ip is the ionization potential, k is the canonical momentum (measured on the de- 
tector), and A(t) = — (F/w) sunp. 



According to Equation (2.41 ), the probability of one-electron ionization F can be written 



as 



r oc exp [— 2ImS' (A;||, A;_L,/p)] 



S {k\\,k±,Ip) 



k) 



(^l[h+A{t)r+h^'+ip^ dt, 



(3.2) 
(3.3) 



where S is the action. Equation (2.42) can be rewritten in terms of S as 



S {kii,k±,Ip) = 0. 



(3.4) 
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Note that the analogy between the saddle point S-matrix calculations |117] . where transi- 
tions are calculated using stationary points of the action, and the adiabatic approach can 
be seen from Equations (3.2) and (3.4). The transition point is given by 

h 



utn = Arcsin 



where 7 is the Keldysh parameter 



7 




(3.5) 



7 



f' 



2U„ 



and Up = {F/2u)'^ is the ponderomotive potential. In order to extract the imaginary and 
real parts of this solution, the following equation |118] can be used 



+ Vi 



Arcsin (x + iy) = 2Ktt + arcsin j3 + i\n 
where K is an integer and 

{^} = lV¥+W+¥±lv^¥^W+ 



(3.6) 



Using Equation (3.6) in Equation (3.2), we obtain 



r(7, k\i, k±) oc exp 
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(3.7) 



where 

/(7,^|hfc±) 



arccosh a — \/ — 1 





2 , a[l 



2/32 



2 

7 



'2L, 



1 



Note that a > 1. It must be stressed here that no assumptions on the momentum of the 
electron have been made. However, Equation (3.7) has an exponential accuracy because 
the influence of the Coulomb field of a nucleus cannot be accounted for by the strong 
field approximation (SEA). The correct exponential prefactor has been obtained within 
the Perelomov-Popov-Terent'ev (nil [I20l [Ml 1122] (PPT) approach. 



Similarly to the Yudin-Ivanov formula |3T], Equation (3.7) is valid if the strength of 
the laser field F depends on time, F — )■ Eog(t), where the envelope g{t) of the pulse is 
assumed to be nearly constant during one-half of a laser cycle. 
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3.2 Connections with Previous Results 



In this section, Equation (3.7) is applied to some special cases in order to establish con- 



nections with previously known results. 

In the case of zero final momentum (fc = 0), we have that a = a/1 + 7^ and /3 = 0. In 
this limit we recover the original Keldysh formula [30\ 



/(7,0,0) 



1 H ) arcsinh7 

27^ 



27 



In the tunneling limit (7 <^ 1) the following formulas can be obtained. Expanding 
the function /(7, /cy, k±) in a Taylor series up to third order with respect to 7 and setting 
k± = 0, we obtain 

r(7,L,,0) ^r(7,0,0)exp 



3u ^ 



Equation ( 3.8[ ) has been derived by a classical approach in Reference |123] (see also Ref- 
erence |114j ). Discussions regarding the physical origin of Equation (3.8) are presented in 
Reference |117j . Performing the same expansion and setting = 0, we obtain 



r(7, 0, k±) oc exp 



2 {ks_^ + 2Ip) 
3F 



3/2' 



(3.9) 



This equation has been derived in Reference |114] . For small values of kj_, Equation (3.9) 
can be approximated by 



r(7,0,fcx) ^r(7,0,0)exp 



/2Ipkj 



(3.10) 



Let us fix fcx = and continue working in the tunneling regime. For the case of high 
kinetic energy, such as k^\'^/2 > 2Up and \ k\\^ / [AUp) — 1 ^ 7, we obtain 



fell 



a 



4f/„ 



/3 



and the ionization rate F is given by 



F oc exp 



2Ur, 



2U„ 



+ 1 arccosh 



4f/„ 



AUp \AUp 



(3.11) 
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Equation (3.11) has been obtained in Reference 



Calculating the asymptotic expansion of the function /(7, /c||,0) for jl ^ 2f/p, we 
obtain 



r oc 



F2 exp(3) 



{k\i^/2 > 2t/p, ki_ = 0). 



(3.12) 



Equation (3.12) has been obtained for tunneling ionization in Reference pi]. Here, we 



have proved that Equation (3.12) is valid for arbitrary 7. A similar formula can be derived 
for k±^Q 



r oc 



exp(l) 



{k^^/2 > 2Up, /c|| = 0), 



4^2^x2 

which is also valid for arbitrary values of the Keldysh parameter 7. 



(3.13) 



Consider the asymptotic expansion of Equation (3.7) for large values of Ip {Ip ^ k^/2). 
In this case a and /3 can be approximated by 



a 



Using these equations, we obtain 

/(7,A;|i,A;x)^/(7,0,0) 



/3 



7 



r 



/2L 



2L 



■arcsinh7 — 



7 



1 + 72 2Ip 



(3.14) 



Equation (3.14) has been reached within the PPT approach (see also Reference [33]). 



Goreslavskii et al. |36| have derived an expression for the spectral-angular distribution 
of single-electron ionization without any assumptions on the momentum of the electron. 
However, they have summed over saddlepoints, i.e., the contribution from previous laser 
cycles has been taken into account. On the contrary, we have not performed any summation 
because we are interested in the most recent contribution to ionization. Therefore, our 
result does account for the phase dependence of the ionization rate, unlike that of Reference 



To make the phase dependance explicit in Equation (3.7), we apply the substitution 

ki\ k\\ - A{t). 



(3.15) 



The analytical expression for the ionization rate as a function of a laser phase when k 



has been achieved by Yudin and Ivanov [31]. Thus, Equation (3.7) is seen to be a 
generalization of the Yudin-Ivanov formula. 
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Note that generally speaking, there is no unique and consistent way of defining the 
instantaneous ionization rates within quantum mechanics, and such a definition is a topic 
of an ongoing discussion (see, e.g.. References |124[ 1125] and references therein). However, 
the instantaneous ionization rates are indeed rigorously defined within the quasiclassical 
approximation (e.g., the Yudin-Ivanov formula), and we have employed this approach. Al- 
ternatively, one can approximate the instantaneous ionization rates by the static ionization 
rates at each point in time using the instantaneous value of the laser field. 



Lastly, we illustrate Equation (3.7) for the case of a hydrogen atom. The single-electron 
ionization spectra in the multiphoton regime (7 ^ 1) and in the tunneling regime (7 <^ 1) 
are plotted in Figure |3.1[ a) and Figure |3.1[ b) respectively. One concludes that the smaller 
7, the more elongated the single-electron spectrum. We can notice that the maxima of 
both spectra are at the origin. Nevertheless, a dip at the origin has been observed experi- 
mentally |126[ I127i I128j in the parallel-momentum distribution for the nobel gases within 
the tunneling regime, and afterwards it has been investigated theoretically in paper |129] 
and references therein. However, such a phenomenon is beyond Equation (3.7). Recently, 



Formula (3.7) was employed to interpreted experimental measurements of perpendicular 
momentum distributions of photoelectrons [130J - substantial, but not total, agreement 
was observed (see Figure 3.4). 



The phase dependence of ionization for different initial momenta, recovered by means 
of Equation (3.15), is illustrated in Figure 3.2 for selected positive momenta. The curves 
for negative momenta are mirror reflections (through the axis = 0) of the corresponding 
positive curves. Figure 3.3 shows that the cutoff of the single-electron spectrum in the 



tunneling regime (the dashed line) corresponds exactly to the kinetic energy 2Up, which 
is the maximum kinetic energy of a classical electron oscillating under the influence of a 
linearly polarized laser fleld. 
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Figure 3.2: (Reprinted from Reference |116] . Copyright (2008) by the American Physical 
Society.) The dependence of the photoelectron spectrum of a hydrogen atom on the instan- 
taneous phase of a laser field: the plot ofT('j,k + ^ sin0, 0) [Equation (3.7)] for a hydrogen 
atom at 800 nm; (a) at 1 x 10^^ W/cm^ (7 = 3.4); (b) at 6 x 10^^ W/cn? (7 = 0.4). 
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Figure 3.3: (Reprinted from Reference Copyright (2008) by the American Physical 

Society.) The longitudinal momentum distribution of photoelectrons of a hydrogen atom: 
the plot of r(7, , 0) /r(7, 0, 0) [Equation 



line: 1 x 10^^ W/cm^ (7 



3.7)] for a hydrogen atom at 800 nm; the solid 
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Figure 3.4: (Reprinted from Reference |130] with kind permission of C. Smeenk. Copyright 
(2010) by the American Physical Society.) (a) Width of the experimentally measured 
momentum distribution in Reference |130j as a function of the strength of the laser field 
for Ar at 800 and 1400 nm. The experimental data are compared with predictions of 

Results obtained by means of these equations are marked as 
, correspondingly, (b) Width of the experimentally measured 



Equations (3.7) and (3.10) 
and 



"non-adiabatic" and "ADK 
momentum distribution as a function of the field for Ne. 
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Chapter 4 

Nonsequential Double Ionization 
below Intensity Threshold: 
Contributions of the 
Electron-Electron and Electron-Ion 
Interactions 



4.1 Formulation of the Problem 

Our model of nonsequential double ionization (NSDI) complements earlier theoretical work 
on calculating correlated two-electron spectra using the strong-field S-matrix approach 
[66] . The key theoretical advance of this work is the ability to include nonperturbatively 
all relevant interactions for both active electrons: with each other, with the ion, and with 
the laser field. Electron-electron and electron-ion interactions are included on an equal 
footing. Our model ignores multiple recollisions and multiple excitations developing over 
several laser cycles, such as those seen in the classical simulations This simplification 
is particularly adequate for the few-cycle laser pulses, as demonstrated in the experiment 
|131j . According to this experiment, multiple recollisions are noticeably suppressed already 
for 12 fs pulses at A = 800 nm. For 6-7 fsec pulses at A = 800 nm, such simplification is 
justified. 



The process of NSDI is shown by the Feynman diagram in Figure 4.1, The system 
begins in the ground state \gg) at time tj. At an instant tf,, intense laser field promotes the 
first electron to the continuum state |k); the second electron remains in the ground state of 
the ion \g^). Recollison at tr frees both electrons. The symmetric diagram where electrons 
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Figure 4.1: (Reprinted from Reference |132j . Copyright (2009) by the American Physical 
Society.) The diagram of NSDI within the SEE mechanism. 



1 and 2 change their roles is not shown, but is included in the calculated spectrum. 



Let us apply the adiabatic approach [Equation (2.41)] to the two-electron process un- 
der consideration. The NSDI within simultaneous electron emission (SEE) has two stages, 
namely ionization of the first electron and the recoUision. Hence, strictly speaking, the adi- 
abatic approximation has to be applied to each of the two stages, since the total amplitude 
of the process is the product of the ionization amplitude and the recollision amplitude. 
However, it is the second (recoUisioin) amplitude that is responsible for the shape of the 
correlated spectra. The first amplitude only gives the overall height of the spectra, as it 
determines the overall probability of the recollision. Since at this stage we are only inter- 
ested in the shape of the correlated spectra, we omit the ionization amplitude from this 
discussion. 

As a zero approximation, we define Ei{ip) and Ef{ip) for the second part of NSDI 
without the Coulomb interaction. Before the recollision at the moment ipr = ootr (Figure 



4.1), one electron is bound and another is free. The classical energy of the system before 



(Pr IS 

where Ip denotes the ionization potential of the ion and A(t) = — (F/w) sin(u;)f:) is the 
vector potential of a linearly polarized laser field. The time of "birth" (ionization) for the 
first electron ipi, = utb{ip) is the standard function of the instant of recollision ipr, which 
is obtained from the saddle-point S-matrix calculations in Appendix |Aj In Equation (4.1) 



we have assumed that the recolliding electron has been born at (pi)((pr) with zero velocity. 
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After the recollision, both the electrons are free and the energy of the system is 



i^/(^) = ^[kl + A((p)]2 + i[k2 + A(^)]^ 



(4.2) 



where ki 2 are the asymptotic kinetic momenta at — ?■ 00 of the first and second electrons, 
correspondingly. 



Now, substituting Equations (4.1) and (4.2) into Equation (2.41), we obtain the corre- 



lated spectrum standard for the strong field approximation (SPA), 



r5irA(ki,k2) OC exp ImS'sirA(ki,k2 



(4.3) 



5'sFA(ki,k2) 



Reip\ 



,0 2 



[ki + A(^)]V-[k2 + A(^)]^ 



-^[A(^)-A($(7;^))]V/f 



where the phase of "birth" (ionization) $(7; (p) corresponding to the recoUison phase cp 
and the transition point ip^ are defined by Equations (A. 11) and (A. 8) in Appendix [A| and 
7 is the Keldysh parameter for the ion [see Equation (A. 4) . 



The major stumbling block is to account for the electron-electron and the electron-ion 
interactions on the same footing, nonperturbatively. To include these crucial corrections, 
we have to include the corresponding Coulomb interactions into Eij{ip). With the nucleus 
located at the origin, the electron-electron and the electron-core interaction energies are 



l/|ri2(t)|, Ki!'') = -2/|ri,2(t)|, 



(4.4) 



correspondingly. Here ri2(t) = ri(t) — r2(t) and ri 2(t) are the trajectories of the two 
electrons. 

However, we immediately see problems. The corrections depend on the specific tra- 
jectory, and one needs to somehow decide what this trajectory should be. Note that the 
classical trajectories ri 2(t) in the presence of the laser field and the Coulomb field of the 
nucleus may even be chaotic. The solution to this problem has already been discussed 
in the Perelomov- Popov- Terent'ev \V20\ \V21\ 1122] (PPT) approach for single-electron 
ionization. In the spirit of the eikonal approximation, these trajectories can be taken in 
the laser field only |12H 1122] I133 l|. so that they correspond to the saddle points of the 
standard SEA analysis. Not surprisingly, in the SEA these trajectories start at the origin, 



ri,2(t)= + A{T)]dr. 

J to 
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(1 2) 

However, here we run into the second problem: the potentials V^e and \4„' are singular. 
Consequently, the integral in Equation (2.41) is divergent and the result is unphysical. 
Therefore, such implementation of the Coulomb corrections requires additional care. 

Appendix [B] describes an approach that deals with these two problems, both defining the 
relevant trajectories and removing the divergences of the integrals. Summarizing the results 
of Appendix |B| we conclude that these problems are overcome by simply using the SFA 
trajectories and soften versions of the Coulomb potentials (4.4). Evidently, there are many 
ways to soften (i.e., to remove the singularity) of the Coulomb potentials; nevertheless, 
final results qualitatively are not affected by such freedom. 



4.2 The Correlated Two-electron Ionization within the 
SFA 



In this section, we analyze the correlated spectrum of the NSDI by using the SFA. In the 
next sections, we will improve the SFA result by employing the perturbation theory in 
action with the SFA result as the zero-order approximation. 



Ignoring the Coulomb corrections in Equation (B.9) and performing the saddle-point 



calculations described in Appendix |A} we reach the usual SFA expression for the correlated 
NSDI spectra - Equation (4.3). 



To illustrate the SFA results, we plot the two-electron correlated spectrum for a system 

In Figure 4.2[ a) we set 7 = 0. Such an SFA spectrum 
—0.78 a.u., which is the maximum of the vector potential 



with Ip of Ar in Figure 



4.2 



has a peak at k\\-^ 

—F/uj ^ —0.78 a.u. The last fact has the following interpretation: NSDI is most efficient 
when the velocity of the incident electron is maximal. This is achieved near the zero of 
the laser field, E(y9) = Fcosy), and the maximum of A((y9). An electron liberated at this 
time could acquire the final drift velocity ~ —F/uo. However, including the correct value 
of the Keldysh parameter 7 not only substantially shifts the peak position [Figure 4.2[ b)], 
but also lowers the maximum by nine orders of magnitude. 



4.3 Electron-Electron Interaction 

In this section, we demonstrate the changes in the correlated spectrum due to the electron- 
electron repulsion. 

Coulomb corrections to the single-electron SFA theory were first introduced by PPT 
using the quasiclassical (imaginary time) method (for reviews, see References [351 1113] ). 
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Figure 4.2: (Reprinted from Reference |132] . Copyright (2009) by the American Physical 



Society.) Correlated two-electron spectra (4.3) of Ar (linear scale) within the SEA at 
7 X 10^3 w/cm2, 800 nm {k^^ = k^^ = 0) (a) 7 = 0; (b) 7 = 1.373. Maxima of figures 
correspond to probability densities: (a) 1.7 x 10^^, (b) 2.9 x 10^^^. 



More recently, further improvements to this method have been considered in References 
|134l I135[ 1136] . These improvements considered not only subbarrier motion in imaginary 
time, but also the effects of the Coulomb potential on the phase of the outgoing wave packet 
in the classically allowed region. These improvements allowed the authors of References 
|134l I135[ I136j to obtain quantitatively accurate results not only for ionization yields, 
but also for the above threshold ionization spectra of direct electrons (i.e., not including 
recollision). An alternative, but conceptually similar, approach is the strong- field eikonal- 
Volkov approach |133j (SF-EVA). Unlike the two previous methods, the SF-EVA allows a 
simple treatment of the electron-electron and electron-ion interaction in the two-electron 
continuum states. 

According to the SF-EVA, the contribution of the interaction potentials is calculated 
along the SEA trajectories. 



1 



ri,2(<^) = - / [ki,2 + A(0)] d(j). 

Note that at the moment of recollision (/?°, the electrons are assumed to be at the origin, 
ri,2(v^r) = 0- However, this does not cause any divergence since according to Equation 
(B.9) we have to use the regularized potential Ke.ing- 



From Equations (B.l) and (B.2), the potential energy of electron-electron repulsion 
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along these trajectories is given by 



Vee,lng{'-P) 



1 — exp 



ruif) = \ri{(f) - r2{(f)\ 



' 12 



(^ill-^b) 



+ 



1 2 



(4.5) 



As discussed in Appendix B , the parameter rjg'' is set to the ionic radius, r^^2 — ^/^p 

The correlated spectrum, which accounts for the electro-electron interaction, has the 
form 



,(0) 



r(2) 



ree(ki, k2) OC CXp ( Im [SsFAlki, k2) + 5'ee(ki, k2) 



(4.6) 



S'ee(ki,k2) 



Vee,lng{<^)d'p. 



Figure 4.3 shows the contribution of electron-electron repulsion to the spectra of NSDI 



for an atom with Ip of Ar (for experimental data see Reference [75] and Figure 4.5). 



Comparing Figures 4^ and |4.3[ we readily notice a dramatic influence of electron- 
electron interaction on the correlated spectra. Electron-electron repulsion splits each SFA 
peak into two peaks because, due to the Coulomb interaction, two electrons cannot occupy 
the same volume. Note that the larger the difference between the perpendicular momenta 
of both the electrons, the closer is the location of the peaks. 



4.4 Electron-Ion Interaction 



Now we include the electron-ion attraction. The potential energy of electron-ion interaction 
for the case of two electrons and a single core, after partitioning (B.l) and (B.2), is 

2 



1 — exp 







^11 1,2' 



+ — (cosv? - cosv3°) ] + [k±i 



(4.7) 



As far as the parameter ro is concerned, we set it equal to r^^2 = ^l^v''' ■ 



r(2) 
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Figure 4.3: (Reprinted from Reference |132] . Copyright (2009) by the American Physical 
Society.) Role of electron-electron interaction. Correlated spec tra of Ar (linear scale) 



at 7 X 10^^ W/cm^, 800 nm are calculated using Equation (|46| with rj"^ = 0.985 a.u. 
(7 = 1.373). Electron-core interaction is not included. Spectra are shown for different 
values of (in a.u.) for both electrons: (a) k^^ = ^±2 = 0; (b) k±i = 0,k±2 = 0.2; 
(c) fcxi = -0.l,k^2 = 0.2; (d) k^, = -0.2,fcx2 = 0.3; (e) k^, = -0.3,k^2 = 0.3; 
(f) k±i = — 0.5, A;_L2 = 0.5. Maxima of figures correspond to probability densities: (a) 



9.1 X 10~i^ (b) 6.5 X 10-^^, (c) 8.36 x 10"^^ (d) 8.4 x 10"^^ (e) 9.1 x 10-^^ (f) 2.8 x 10' 
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Finally, the correlated spectrum of NSDI, which takes into account both the electron- 
electron and electron-ion interactions, reads 



S'e„(ki,k2) 



ree+en(kl, k2) OC exp ( lui [S s F A^^l , '^^2) + Seei^l, ^^2) + 5'e„(ki, k2)] ) , (4.8 



ReipO 



To illustrate the influence of the electron-ion attraction, we have plotted the correlated 



spectra of Ar in Figure 4.4 for different perpendicular momenta. Comparing Figures 4.3 and 



4.4[ we conclude that the larger the difference between the perpendicular momenta of the 
two electrons, the larger is the contribution of the electron-ion interaction. Furthermore, 
accounting for electron-ion attraction increases the probability of NSDI by 15 orders of 
magnitude. This occurs because, as in the case of single-electron ionization, electron- 
core interaction significantly lowers an effective potential barrier. We can also conclude 
that correlated spectra pictured in Figures 4.4 c), 4.4[ d), and 4.4[ e) have the biggest 



contribution to the total probability of NSDI, which is an integral of the probability density 
over momenta of both the electrons. Note that, on the one hand, the maximum of the 
probability density shown in Figure |4.4Kf) is the largest among those presented in Figure 



4.4[ and on the other hand, this maximum is localized in a few pixels; therefore, the 
integral contribution of Figure [4!4| (f) to the total probability is smaller than Figure 4.4[ e). 
Additionally, as one would expect, further increasing of k± leads to a decrease of probability 



density. The correlated spectra in Figure 4.4 agree with the experimental data [7H| [7B] in 



quadrants one and three (see Figures 4.5 and 4.6). The considered diagram (Figure 4.1) 



does not contribute to signals in quadrants two and four. Note that taking into account a 
nonzero value of 7 is vital to achieve agreement with the experimental data. 



From Equations (4.5) and (4.7), we can notice that if ~^ ^ ^^"^ '"o ~^ 00, the 
Coulomb corrections Vexing and Ven,ing vanish, and the SFA result is recovered. Therefore, 
we conclude that the radii r^y} and tq contain the information about the initial position 
of electrons after they emerged in the continuum. Obviously, the intra-electron distance 
should be on the order of an ion radius. 



4.5 Conclusions 

The analytical quantum-mechanical theory of NSDI within the deeply quantum regime, 
when the energy of the active electron driven by the laser field is insufficient to coUisionally 
ionize the parent ion, has been formulated based on the adiabatic approach. On the whole, 
the presented model qualitatively agrees with available experimental data [721 [751 US] • We 
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Figure 4.4: (Reprinted from Reference |132] . Copyright (2009) by the American Physical 
Society.) Role of the electron-core and electron-electron interactions. Correlated spectra 
of Ar (linear scale) at 7 x 10^^ W/cm^, 800 nm are calculated using Equation (4.8) with 



rf^ = fo = 0.985 a.u. (7 = 1.373). Both electron-elecron and electron-core interactions 
are included. Spectra are shown for different values of k± (in a.u.) for both electrons: (a) 
k±i = k^2 = 0; (b) fcxi = 0,fcx2 = 0.2; (c) k±i = -0.1,^^2 = 0.2; (d) fc^i = -0.2, A;x2 = 
0.3; (e) k±i = — 0.3, A;_l2 = 0.3; (f) k±i = —0.5,k±2 = 0.5. Maxima of figures correspond 



to probability densities: (a) 
4.4 X 10-^ (f) 7.7 X 10"^ 



2.8 X 10-^ (b) 3.3 X 10-^ (c) 6.1 x 10-^ (d) 1.5 x 10-^ (e) 
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Figure 4.5: (Reprinted from Reference [75j with kind permission of Y. Liu. Copyright 
(2008) by the American Physical Society.) Experimentally measured correlated spectra of 
Ar at 800 nm (a) 9 x lO^^ W/cm^ [72], (b) 7 x lO^^ W/cm^ [73], and (c) 4 x 10^^ W/cm^ 
[75]. 




Figure 4.6: (Reprinted from Reference [76j with kind permission of Y. Liu. Copyright 
(2010) by the American Physical Society.) Experimentally measured correlated spectra of 
(a) Ar at 3 X 10^^ W/cm^, 800 nm and (b) Ne at 1.5 x 10^^ W/cm^ 800 nm in Reference 

m- 
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have defined the quantum-mechanical phase of birth of the active electron (A.ll), which 
accurately accounts for tunnelling of the recoUiding electron in the regime where both the 
phases and ipi, are complex. Moreover, it has been demonstrated that ignoring such a 
contribution of tunnelling of the active electron fails to agree with the experimental data. 

Furthermore, our results show that any attempt to interpret NSDI spectra in this regime 
in terms of a simple SFA-based streaking model would lead to wrong conclusions on the 
relative dynamics of the two electrons. 

The contributions of the electron-electron and electron-ion interactions have been an- 
alyzed. Both play an important and distinct role in forming the shape of the correlated 
spectra. 

The presented model is not able to reproduce the correlated spectra obtained exper- 



imentally [75l 176] in quadrants two and four (see Figures 4.5 and 4.6). It is because 
the considered process of SEE, when two electrons detach simultaneously from the atom, 
does not contribute to that area. However, it is widely accepted that the anticorrela- 
tion of the electrons in those parts of the correlated spectra is formed due to recoUision- 
induced excitation of the ion plus subsequent field ionization of the second electron (RESI) 
|137l I138[ [THl [THl EH [82] , and it should be noted that this mechanism has been also observed 
in classical simulations [581 EDI [83] . 
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Chapter 5 

Nonsequential Double Ionization 
below Intensity Threshold: 
Anticorrelation of Electrons without 
Excitation of the Parent Ion 



Is excitation of the parent ion indeed necessary to explain the electron anticorrelation in 
nonsequential double ionization (NSDI) below intensity threshold (BIT)? 

We show that this is not always the case. Our conclusions are based, first, on model 
ab-initio calculations showing that the anticorrelation of the electrons exists even if the 
ion has only a single bound state. Second, we present a simple analytical model based 
on the assumption that both the electrons are ejected simultaneously (the time of return 
of the first electron coincides with the time of liberation of the second electron, i.e., the 
ion has no time to be excited). An advantage of this model is that it allows for a simple 
analytical solution in closed form. In a certain range of parameters, the correlated two- 
electron spectrum obtained within this model exhibits the anticorrelation of the electrons. 
This novel mechanism of simultaneous electron emission can produce the anticorrelation 
of the electrons in the deep BIT regime. 



5.1 An Ab Inito Evidence of the New Mechanism 



Consider the model Hamiltonian for a system of two one-dimensional electrons 



^^^^ = -2(a^ + a^)+^(^0 + ^(^2)-n^2-a;i) + (a;i + a;2)F4t), (5.1) 
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where V{x) = — 4exp(— 3x^) is the prototype for the potential of the electron-core at- 
traction, —V{x2 — Xi) is the prototype of the electron-electron repulsion, and F^it) = 
Ff{t)sm{ojt) {F = 0.05 and u = 0.6) - the laser pulse, where f(t) is a trapezoid with 
one-cycle turn-on, six-cycle full strength, and one-cycle turn-off. The potential V{x) is 
chosen such that the one-particle Hamiltonian, Hion = —d'^/{2dx^) + V{x), supports only 

(2) 

one bound state with the ionization potential Ip = 2.11. 

To make analysis more transparent, we substitute the original problem (NSDI) by the 
corresponding problem of laser-assisted scattering, i.e., we simply discard the first step 
of NSDI - liberation of the first electron. In other words, instead of assuming that the 
two-electron system initially is in the ground state, we assume that the first electron is 
an incident wave packet and the second electron is in the single bound state of the ion. 
This modification is very useful, as it allows us to eliminate all other possible interactions 
and processes except the three major components - the electron-electron repulsion, the 
electron-core attraction, and the electron-laser field interaction. 

We solve numerically, by means of the split-operator method, the time-dependent 
Schrodinger equation, 

id<i/{xi,X2;t)/dt = ^(t)^(xi, xg; t). (5.2) 
The initial unsymmetrized wave function takes the form 

X2; 0) = exp [-{xi - /i) V(2fT^)^] ^ix2), (5.3) 

where ax = 2, fi = 7(Tx, and Hianip{x) = —Ip tp{x). The parameters are selected such that, 
Ip^^ — {F/uY /2 ^ 4^, i.e., the bound electron needs to absorb at least four photons to be 
liberated. The wave function then is symmetrized, i.e., we assume that the spins of the 
electrons are antiparallel. 

The wave function of the model system in the coordinate and momentum representa- 



tions is pictured in Figure |5.1[ The wave function in the absence of the laser produces 
only maxima on the axes xi^2 corresponding to one electron bound and the other one free; 
hence, this part of the wave function, which also shows up when the laser file is turned on. 



should be ignored because it does not correspond to NSDI. From Figure |5.1[ we observe 
that the two electrons "prefer" to be anticorrelated rather than correlated even though we 
effectively soften the repulsion between them by selecting antiparallel spins. 



5.2 An Analytical Quasiclassical Expression of Corre- 
lated Spectra for the SEE Mechanism 

We shall calculate the correlated spectrum within the adiabatic approximation. 
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As in the previous section, we replace the problem of NSDI BIT by the problem of 
laser-assisted scattering of an electron by an ion. We study the simultaneous electron 
emission (SEE) process, when the moment of collision of the incident electron coincides 
with the moment of ionization of the ion. 

To include electron-electron repulsion, we follow Chapter |4] and apply the standard 
exponential perturbation theory - strong-field eikonal-Volkov approach |133] (SF-EVA). 
Namely, first we find electron energies and trajectories without the electron-electron repul- 
sion. Then, we correct electron action and energies by adding the effect of the electron- 
electron repulsion, which are calculated along the zero-order trajectories. 

Thus, in zero order, we define Eij{ip) without the electron-electron repulsion. Since 
before ionization, one electron is free and the other is bound, the total classical energy of 
the system before the collision is 

E^{^) = ^[P + AM]^-/f, (5.4) 

where p is the canonical momentum (i.e., the kinetic momenta aX (f = ±00) of the incident 
electron, Ip is the ionization potential of the ion, and A{ip) = — (F/w) simp is the vector 
potential of a linearly polarized laser field. After collision both the electrons are free, and 
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the classical energy of the system reads 

1 



Efi^) = -[kl + A(v.)]V-[k2 + A(v.)]^ 



(5.5) 



where ki_2 are canonical momenta of the first and second electrons. 

Such two electron process is formally equivalent to single-electron strong field ionization 
of a quasiatom (within a pre-exponential accuracy). This statement is manifested by the 
following identity 



1 



1 

UJ 



(5.6) 



where the right hand side of Equation ( 5.6[ ) is the action of single-electron ionization 
within the strong field approximation (SFA) [Equation (3.3)], ,Jt = fcy-^ + ^Il2 ~ 



is the 



effective longitudinal momentum, and 



+ (kj + k2 - p2 - jr2) /2 is the effective 



ionization potential of the quasiatom. Therefore, derivation of the correlated spectra for 
SEE is reduced to calculation of the momentum distribution of photoelectrons after single- 
electron ionization without any restrictions on Jf. The most suitable solution of the last 
problem for our current discussion is given by Equation (3.7). However, this equation is 
valid only for positive J^^; hence, the case of J^p < being of interest for SEE, must also 
be considered. 



Substituting Equations (5.4) and (5.5) into Equation (2.41) and taking into account 

(5.7) 



the previous comments, we obtain 

r(ki, ka) oc exp (-2 G/u + Ve 



where 



G 

y± 

9 



arcsinx_ -|- iarccoshx_|. if J^p < 
arcsin ?/_ + iarccosh ?/+ if j^p > 0, 

1 ) arccosha;+ — ^^x\ — 1 

j arccosh — ^y\ — \ 

1^7(77 -l) + l|/2±|(7(r7-l)-l|/2. 



2r\X— 
9 



+ X, 



1-2x1 



if J^p < 0, 
if ^p > 0, 



(5.8) 
(5.9) 



2\.yp\/F, ri = je/^2\yp\. 



In Equation (5.7), V^e denotes a crucial correction due to the electron-electron repulsion. 



-— Im 
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-1X1 |ki 



Rc 1^0 

k2 



\vx{ip) -V2{^)\ i™|ki-k 



rim 



(5.10) 
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Figure 5.2: Correlated two-electron spectra of Ar (linear scale) given by Equation (5.7) at 
800 nm and 1 x 10^^ W/cm^ (p^ = = ^±2 = 0) for different momenta of the incident 
electron: (a) py = 0.1, (b) py = 0.2, (c) py = 0.25, (d) py = 0.4. 
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Figure 5.3: Correlated two-electron spectra of Ar (linear scale) given by Equation (5.7) at 
800 nm (p = and = k±2 = 0) for different intensities of the laser field: (a) 3 x 10^^ 
W/cm2, (b) 2.7 X 10^3 w/cm^, (c) 2.5 x lO^^ w/cm^, (d) 2.3 x lO^^ w/cm^, (e) 2 x lO^^ 
W/cm2, and (f) 1 x lO^^ w/cm^. 
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The model of NSDI BIT presented here [Equation (5.7)] is similar to the one developed 
in Chapter |4| However, there is an important difference. In the correlated spectra (5.7), 
the momentum of the incident electron p is a free parameter, whereas the correlated 
spectrum given by Equation (4.8) does not have this freedom - the canonical momentum 



of the recoUiding electron is fixed as a function of the phase of recoUision [see Eq. (A. 11)] 



Therefore, Equation (5.7) allows us to establish the link between values of the canonical 



momentum of the rescattered (incident) electron and specific portions of the correlated 



spectra. The dependance of the correlated spectra (5.2) on the momentum of the incident 



electron is presented in Figure 5.2 



However, more interesting question is the dependence of the correlated spectrum (5.7) 



on the intensity of the laser field, which is pictured in Figure 5.3 The parameters used 



to plot Figure 5.3[ a) coincide with the parameters employed in the recent experiment 
[76] (see also Figure 4.6). Yet, the experimentally measured correlated spectrum exhibits 

The reason 



peaks in the second and fourth quadrants, which contradicts Figure |5.3K a) 
of such a disagreement is that the anticorrelation in these experimental data is due to 
the recoUision-induced excitation of the ion plus subsequent field ionization of the second 
electron (RESI) mechanism. As the intensity lowers, the peaks in the correlated spectra 
of SEE shift to the second and fourth quadrants [see Figure 5.3[ f)], i.e., the SEE process 
leads to the anticorrelation of the electrons in the deep BIT regime. This can be explain 
intuitively in the following way: The lower the intensity, the smaller the canonical mo- 
mentum of the returning (incident) electron. Since the canonical momentum of the system 
is approximately conserved, we obtain ~ ki + hence, ki ^ — k2. Having absorbed 
a necessary number of photons, both the electrons emerge in the continuum where they 
experience strong electron-electron repulsion that pushes them apart. Therefore, during 
SEE, the electrons gain momenta because of the electron-electron repulsion. Indeed, if the 
electron-electron interaction is "turned off" in Equation (5.7) by setting Vee = 0, then we 



obtain a single peak centred at the origin instead of the two peaks visible in Figure |K3](f). 
Recall that contrary to SEE, the electrons gain kinetic energy by oscillating in the laser 
field in the course of RESI. 



5.3 Conclusions 

We have demonstrated that the mechanism of simultaneous electron emission, when the 
time of return of the rescattered electron is equal to the time of liberation of the bounded 
electron, can be responsible for the anticorrelation of the electrons during NSDI in the deep 
BIT regime [see Figure [5l3| (f)]. The SEE process significantly differs from RESI because 
it does not require the excitation of the ion to explain the anticorrelation of the electrons 
observed in the two-electron correlated spectra. SEE and RESI are by no means mutually 
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exclusive; they both contribute. It will be interesting to study quantitatively the relative 
contribution of SEE and RESI to NSDI in the deep BIT regime. 
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Chapter 6 



Shapes of Leading Tunnelling 
Trajectories for Single-electron 
Molecular Ionization 

If the frequency of the laser field is very low, then it is a good approximation to reduce 
the original time dependent problem to the time independent one where the laser field 
is substituted by a static field. In such a picture, ionization processes are realized by 
quantum tunnelling. From the interpretational point of view, it is advantageous to use the 
language of quasiclassical trajectories to calculate tunnelling rates. Since the quasiclassical 
approximation in the original form cannot be readily utilized in the three-dimensional 
space, the assumption that all relevant quasiclassical trajectories are close to linear is very 
important, because it reduces the original three-dimensional problem to the effective one- 
dimensional problem. Therefore, it is not surprising that such an approximation became 
almost a tacit assumption in strong filed physics, especially in the case of single-electron 
molecular ionization. In fact, we have already used a very similar assumption in the time 
dependent cases when we calculated the Coulomb corrections - The trajectories of electrons 
in the combined field of the core and the laser have been assumed to be merely strong field 
approximation (SFA) trajectories perturbatively corrected by the Coulomb field. 

How good is the hypothesis of linearity of tunnelling trajectories in strong fields? To 
answer this question, we first need to pose it clearly. 

6.1 Mathematical Background 

The instanton approach is one of the methods for description of tunnelling |139j . It can be 
introduced as a result of application of the saddle point approximation to the modification 
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of the Feynman integral obtained after performing the transformation of time t — > —ir 
to "imaginary time" r (i.e., the Wick rotation). This technique has turned out to be 
tremendously fruitful in many branches of physics and chemistry (see, e.g.. References 

We shall reiterate the main steps in deriving the instanton approach. Let us consider 
a quantum system with the Hamiltonian 



H = -A/(2m) + [/(x) 



(6.1) 



where A is the n-dimensional Laplacian and x is an n-dimensional vector. The Feynman 
integral representation of the propagator reads |144j 



5[x] = 



|x,) = N I iF[x(t)]e*^["W] 

-to 



(6.2) 



J^(x,x)(it, if(x,x) 







2m 



t/(x) 



where the path integral sums up all the paths that obey boundary conditions x(0) = Xj 



and x(to) = x/, and x(t) = (h(.{t)/dt. After performing the Wick rotation. Equation (6.2) 
becomes 



S[x(r)] 



(6.3) 



^[x] 



TO 



1 /c/x(rl 



2m \ dr 



+ f/(x(r)) 



where Tq = Uq and S is called the Euclidian action. Hence, one can say that the transition 



from Equation (6.2) to Equation (6.3) is achieved by the following formal substitutions 



t — )■ — zr, x(t) — 7- x(r), x(t) — )■ id'x.{T)/dT. 



(6.4) 



Comparing the actions S and S, one concludes that the motion in imaginary time is 
equivalent to the motion in the inverted potential. In other words, the actions S and S 
are connected by the substitution 



U -U, {E -E). 



(6.5) 



The final step in the instanton approach is the application of the saddle point approxima- 



tion to the Euclidian Feynman integral in Equation (6.3) assuming that tq — ?■ oo 



However, there is a long ongoing discussion |145t I146[ I147[ 1148] whether the instanton 
approach agrees with the quasiclassical approximation for tunnelling; some observations 
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have been made that these two methods may disagree up to a pre-exponential factor. 
Furthermore, as it has been pointed out in Reference |141] . the instanton approach in 



the formulation presented so far [substitutions (6.4)] not only looks like a "highly dubious 
manoeuvre," but also gives no prescription for getting a correct pre-exponential factor. 
Consequently, a natural question aries how this method can be safely used and what the 



meaning of substitutions (6.4) and (6.5) is. 

The mathematical physics community has reinterpreted the instanton approach rigor- 
ously (see, e.g.. References UMl ESI M [ffl [M HSl [M] and references therein), 
and the corresponding rigorous analysis answers both questions. Moreover, this rigorous 
interpretation is extremely useful because it can be implemented as an effective numerical 
method, which will lead to a clear physical picture applicable to a broad class of problems. 
We shall review briefly the cited above works since on the one hand, they are unknown for 
physicists, and on the other hand, they may be challenging to read for non-specialists in 
mathematical physics. 

Historically, the first problem considered within such a framework was "how fast does 
a bound state decay at infinity?" |149[ 1150] (see also Section 3 of Reference |106j ). Let 



us clearly pose the question. Consider the Hamiltonian (6.1) as a self-adjoint operator 
on ^2{^"') ^ the space of square-integrable functions. A bound state ip G ^2(1^") is a 
normalizable eigenfunction of such a Hamiltonian, Hip = Eip. Since the normalization 
integral converges, the bound state ip = "^(x) must vanish as ||x|| — i- 00. Therefore, we 
want to determine how this decay is affected by the potential V. This question can be 
answered very elegantly if we confine ourself to an upper bound on the rate of decay. 

To obtain this upper bound, we need to introduce first some geometrical notions. Let 
M be a real n-dimensional manifold (intuitively, M is some n-dimensional surface). The 
tangent space at a point x G M, denoted by Tx(M), is a real linear vector space M" that 
intuitively contains all the possible "directions" in which one can tangentially pass through 
X. A metric is an assignment of an inner (scalar) product to Tx(M) for every x G M. 

Let X G M" and ^, ^7 G Ty^{M). We define a (degenerate) metric by 

(^,77), = 2m([/(x)-i5;) + (^,77), (6.6) 

where {$,,11) = ■ rj = ^ir/i + . . . + C,nVn is the Euclidean inner product and /(x)+ = 
max{/(x),0}. Following the convention used in mathematical literature, we shall call 



metric (6.6) as the Agmon metric. 

Having introduced the metric, we can equip the manifold M with many geometrical 
notions such as distance, angle, volume, etc. The length of a different iable path 7 : [0, 1] — )■ 
in the Agmon metric is defined by 

LAh) = f m)hif)dt = f[U{i{t)) - E]fm)\\dt, (6.7) 
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where ||^|| = a/ ^) is the Euchdian (norm) length, and ||^||x = a/ ^)x- The path of a 
minimal length is called a geodesic. Finally, t/ie Agmon distance between points x, y G M", 
denoted by Pe{'^, y), is the length of the shortest geodesic in the Agmon metric connecting 
X to y. 

Before going further, we would like to clarify the physical meaning of the Agmon metric. 
Let us recall the Jacobi theorem from classical mechanics (see page 150 of Reference [153j 
and page 247 of Reference |154] ): The classical trajectories of the system with the potential 
U{x.) and a total energy E are geodesies in the Jacobi metric 

((^,r7))x = 2m(E-f/(x))+(^,r7), (6.8) 

on the set {x G ]R"|f/(x) ^ E} - the classical allowed region. The Agmon metric [Equation 



(6.6)] and the Jacobi metric [Equation (6.8)] are indeed connected through the substitution 
(6.5 ). By virtue of this analogy, we conclude that the Agmon distance has to satisfy a time- 
independent Hamilton- Jacobi equation, also known as an eikonal equation, 

|VxPi^(x,y)|2 = 2m(t/(x) - E)+, (6.9) 

where Vx/(x) = {df /dxi, . . . , df /dxn)- In fact, the Agmon distance is the Euclidean ver- 
sion of the reduced action [here, the adjective "Euclidian" means the same as in Equation 



(6.3)]. In other words, the Agmon distance is the action of an instanton. 

Now we are in position to recall upper bounds on a bound eigenstate of the Hamiltonian 



(6.1). First, under very mild assumptions on U (merely, continuity, compactness of the 
classically allowed region, and absence of tunnelling, i.e., the spectrum of the Hamiltonian 
being only real), it has been proven |150] that for an arbitrary small e > 0, there exists a 
constant < < oo, such that 



e2(i-^)''^W|^(x)|Vx^c„ (6.10) 



where p£;(x) =pe(x, 0). Roughly speaking, result (6.10) means that ?/'(x) = O (e"*^^"^^''^'^^^) . 
However, this result can be improved. For any small e > 0, there exists a constant 
< Ce < oo, such that the following inequality is valid under additional conditions of 
regularity of the potential U 

|^(x)| ^ c,e-(^-^^''^W. (6.11) 



Analyzing Equation (6.10) and Equation (6.11), we conclude that the Agmon distance 



from the origin describes the exponential factor of the wave function. Further information 
can be found in References |150[ \152\ 1106] and references therein. We note that lower 
bounds on ground states can also be obtained by utilizing the Agmon approach |149j . 
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We illustrate the power and utility of upper bound (6.11) by deriving upper bounds 
for matrix elements and transition amplitudes in Appendix [Cj The former result is an 
estimate of the modulo square of the matrix element {ipp\ V {ipq) , where ipp and ipq are 
bound eigenstates of the Hamiltonian (6.1) that correspond to eigenvalues Ep and Eg. It 
is demonstrated in Appendix [C] that for an arbitrary small e > 0, there exists a constant 
< Ce < oo, such that 

\{ijp\V\^g)f ^c,j y2(x)e-2(i-^)KW+''^,W]d-x, (6.12) 

which could be interpreted as, 

\{iJp\V\i:g)f = o(^j y2(x)e-2(i-^)hpW+''^.W]d"x^ . (6.13) 



Simplicity of the derivation of Equation (6.13) does not imply its insignificance. On 



the contrary, Equation (6.13) is a multidimensional generalization of the Landau method 
of calculating quasiclassical matrix elements |3] (see also page 185 of Reference |110j and 
References |155[ 1112] ). To the best of my knowledge, such a generalization has not been 
reported before. To prove the one-dimensional version of the Landau method using an- 
alytical techniques (as it is usually done), one deals with the Stokes phenomenon (see, 
e.g.. Reference |156] ): thus, the generalization to the multidimensional case without too 
restrictive assumptions is not obvious. The Agmon upper bounds lead not only to quite a 
trivial derivation, but also to an intuitive physical and geometrical picture. 

Now we explain briefly how these geometrical ideas are generalized to the problem of 
tunnelling (an interested reader should consult References |103l 11041 1105] 1106] and refer- 
ences therein for details and further developments). Let E be an energy of a tunnelling 
particle. We denote the boundary of the classically forbidden region by Se- It is assumed 
that Se consists of two disjoint pieces and (i.e., Se = S^U S^ and S'g n S'g = 0) - 
inside and outside turning surfaces, which are merely multidimensional analogs of turning 
points. Having introduced the concept of the Agmon distance, we naturally introduce two 
related notions: First, the Agmon distance from the surface S'^ to a point x, p£;(x, 5*^), 
as the minimal Agmon distance between the point x and an arbitrary point y G 5*^ [more 
rigorously, p£;(x, S^) = infyg^- p^(x, y)]; second, the Agmon distance between the turning 
surfaces S'^ and S'^, pe{S^, S^) , as the minimal Agmon distance between arbitrary two 
points x e S"^ and y G S"^ [ Pe{Se, S^) = inf^g_5+ pij(x, S^)]. 

In a nutshell, and thus a bit abusing the formulation of the original result |105] , we say 
that for an arbitrary small e > 0, there exists a constant c > 0, such that the tunnelling 
rate, F, (viz., the width of a resonance) in the quasiclassical limit {h — )■ 0) obeys 

T ^ cexp[-2(3E{pE - e)], (6.14) 
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where < p^; < oo and PePe being the leading asymptotic of pe{S^,S~^) when /i — )• 0. 



However, the following interpretation of upper bound (6.14) is sufficient for our further 
applications 

r = O (e-2''-(^^'^^)) , (6.15) 

i.e., twice the Agmon distance between the turning surfaces gives the leading exponential 
factor of the tunnelling rate within the quasiclassical approximation. 



Equation (6.15) is not only of analytical interest, but also is a starting point of an 



efficient numerical method for estimating tunnelling probabilities. The Agmon distance 



between two points, p£(x, y), can be computed by solving numerically Equation (6.9) with 
the boundary condition 

Pi?(y,y) = o (6.16) 

by means of the fast marching method |157l 1158] I159[ [160] . Moreover, having computed the 
solution, one can readily extract the minimal geodesic from a given initial point x by back 
propagating along p£;(x, y), where y is regarded as a fixed parameter; more explicitly, the 
minimal geodesic, g = g(t), is obtained as the solution of the following Cauchy problem 

g = -V4Pij(^,y), g(0)=x. (6.17) 
Such a geodesic can be interpreted as a "tunnelling trajectory." 



A brief remark on types of the solutions of Equation (6.9) ought to be made. Generally 
speaking, an eikonal equation admits a local solution under reasonable assumptions, but 
a global solution is not possible in a general case owing to the possibility of development 
of caustics (see, e.g.. Reference |161j ). Nonetheless, when we talk about a solution of 



Equation (6.9), we actually refer to a viscosity solution because not only it is a global 
solution, but also it has the meaning of distance |159[ 1160] which we originally assigned to 
the function pE- 

In fact, the fast marching method is an "upwind" finite difference method that efficiently 
computes the viscosity solution of an eikonal equation. Note, hence, that the fast marching 
method as well as the other ideas presented and developed in the current work cannot be 
employed to study the influence of chaotic tunnelling trajectories (see Reference |162j and 
references therein). Some implementations of the fast marching method as well as the 
minimal geodesic tracing can be downloaded from References |163l 11641 1165] . 

The Agmon distance from the surface to a point, pEi^,S^), must satisfy Equation 



(6.9). Indeed, p£;(x, S^) is the solution of the boundary problem 

|VxPi?(x,^^)|2 = 2m([/(x)-E)+, p^(y,^-) = 0, VyG5^, (6.18) 
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which can be solved by the fast marching method as welL Finally, the Agmon distance 
between the turning surfaces is computed as min^g^+ pEi^, 5'^) after solving Equation 



( I6.18D . 

The points h E and e G 5*^ such that 

PEiS^,S+)=pE{h,e), (6.19) 

are of physical importance because they represent the points where the particle "begins" 
its motion under the barrier (b) and "emerges" from the barrier (e), correspondingly. 



Moreover, the minimal geodesic (6.17) that connects these points (g(0) = b and g(l) = e) 
is a tunnelling trajectory which gives the largest tunnelling rates - the leading tunnelling 
trajectory. Note, however, that these points as well as the trajectories may not be unique 
in a general case. 

It is noteworthy that a power of the fast marching method in applications to tunnelling 
has already been recognized in chemistry within the context of the reaction path theory 
|166l 11671 11681 1169j . Similarly to the current work, the main object of interest of those 
studies is the reaction path, which is the leading tunnelling trajectory in our terminology. 
Nevertheless, the motivation for the usage of the fast marching method, presented in 
References [166^ 1167^ 1168^ 1169] , is tremendously different from our geometrical point of 
view. 



6.2 Main Results 

In this section, we shall follow a two step program. First, we consider tunnelling in mul- 
tiple finite range potentials, where we prove that leading tunnelling trajectories are linear 
(Theorem [T]). Then, we reduce the case of multiple long range potentials to the previous 
one by employing the fact that a singular long range potential can be represented as a 



sum of a singular short range potential and a continuous long range tail [Equation (6.36) 



Such a reduction allows us to prove that the leading tunnelling trajectories are "almost'' 



linear (Theorem |2]). We note that partitioning (6.36) was put forth by Perelomov-Popov- 
Terent'ev [ng IM II2I1 [122] (PPT), and it is widely used for obtaining the Coulomb 
corrections in strong filed ionization (see References [35l 11131 11351 11341 11251 1136] and ref- 
erences therein). 

Let us introduce some notations. Hereinafter, the dimension of the space is assumed 
to be n ^ 2. The interaction of an electron with a static electric field of the strength 
F is of the form Fxn {F > 0). dA denotes the boundary of the region A. The map, 
miria;^ : M" D A — i- M", selects a point x = min^.^A G A that has the smallest x„ 
component among all the other points from A, assuming that A has such a unique point. 
The projection Px of the point x = (xi, X2, ■ ■ ■ , Xn) is defined as Px = (xi, . . . , E/F). 
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Theorem 1. We study single electron tunnelling (—00 < E <0, F > 0) in the potential 

K 

f/(x) = ^r,(||x-R,||) + Fx.. (6.20) 

Let us assume that 

1. Vj : {0,Rj) (-00,0) and Vj : {Rj,+oo) {0}, Rj > 0, j = 1, . . . , K , are 
differentiable on {0,Rj) and strictly increasing functions, such that Vj{0) = —00 and 
Vj may have a jump discontinuity at the point Rj . 

2. supp V^- = {x G M" I V,(||x - RjII) 7^ 0} is the support of the potential Vj(||x - Rj||), 
supp Vk n supp Vj = 0, Vfc 7^ j and supp Vj fl {x G M" | x„ ^ E/F} = ij), j = 1, . . . , K . 

3. Introduce qj = min^^ 9supp V^- , pj = min^.^ S'^(j), Sj^{j) is defined in Equation 



(6.22). If there exists N, such that 

\\PN - PRtvII < ||q, - PRill, Vj ^ N, (6.21) 

Then, the leading tunnelling trajectory is unique and linear, and it starts at the point 
and ends at PUn, Pe{S'^, S^) = pe{pn, P^n)- 

Proof. The boundary of the classicaUy forbidden region is defined by the equation U (x) = 
E. Consider two cases: 

First, if X ^ UjLi supp Vj then according to assumption 2, the above equation simply 
reads Fxn = E, and thus its solution defines the outer turning surface 

5+ = {xGM"|x„ = E/F}. 

One can see now that the projector operator P projects a point onto S^. 

Second, if x G supp V^- and Vj is continuous at the point Rj, then the equation reads 
Vj(||x — Rjll) + Fxn = E. To proof that the set 

SeU) = {x G supp Vj I V,(||x - R,||) + Fxn = E} (6.22) 

is not empty, we construct the function /j(x) = V,(||x — Rj||) + Fxn — E. Since fj(Rj) = 
—00, we can find a set Aj C suppV^- located close to Rj, such that /j(x) < for all 
X G Aj] correspondingly, since according to assumption 2, x„ > E/F, there exists the 
set Bj C supp Vj of points close to the boundary of supp Vj for which fj is positive. 
In fact, Aj and Bj can be constructed such that ||x — Rj|| < ||y — Rj||, Vx G Aj and 
Vy G Bj. Therefore, the intermediate value theorem guarantees that Sj^{j) 7^ and it 
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"lies between" Aj and Bj. Furthermore, the inner turning surface is = [jf=i S^{j), 
and S'^{j) n S^{k) = 0, Vj k. (Note that the strict monotonicity of Vj{x) assures that 
the set S'^{j) is connected.) Whence, 

Pe{Se, S^) = min {pi?(^^(j), S^)} . (6.23) 



Equation (6.23) means the reduction of the many centre case to the singe centre case under 
the assumptions made. Needles to mention that such a reduction tremendously simplifies 
the analysis. 

The same conclusions are valid if the jump of the function Vj at Rj is not too large, 
so that the equation V^ (||x — Rj||) + Fxn = E has solutions for x G supp Vj. However, if 
the jump is too large, i.e., this equation does not have solutions from the support of the 
potential, then it is natural to set S'^{j) = dsuppVj. 

Consider the single centre case - single electron tunnelling in the potential Uj{x) = 
V^(||x — Rjll) + Fxn- We shall show that this potential is axially symmetric. If x = 
(xi, . . . , Xn), then we introduce IIx = (xi, . . . , x„_i). We can then symbolically write 
X = (IIx, x„). Using this new notation, we obtain 

f/,(x) = V, (^||nx-nR,P+(x„-[R,]j'^ + Fx„, (6.24) 

where [a] denotes the n*'^ component of the vector a. It is readily seen from Equation 



(6.24) that the potential ?7j(x) is invariant under transformations that do not change Xn 
and arbitrary (n — l) dimensional (proper and improper) rotations of the vector IIx about 
the point IIRj . The only invariant subspace of M" under such transformations is the line 

{{URj,Xn)\XneR}. 

Since both regions Sj^{j) and are shape invariant under the axial symmetry trans- 
formations, we may expect that the shortest geodesic connecting these regions ought to be 
shape invariant as well. Thus, one readily concludes that the leading tunnelling trajectory 
should be linear and should connect the points and -PRj 

PE{Pj,PRj) = Pi^(S^(j),5+), (6.25) 

since no other geodesic that connects S^{j) and is shape invariant with respect to the 
axial symmetry transformations. Below we shall present a formal version of this derivation. 

Foremost, we demonstrate that the operation minx„ is defined on the set S^{j), viz., 
that there is a unique point of Sj{j) that has the smallest component x„. Employing the 
method of Lagrange multipliers and taking into account the symmetry of the potential, we 
construct the function 



^l(x„,C, A) = Xn + X 



(6.26) 
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The condition d^i/dc = leads to c = 0. Therefore, pj = min^.^S'£;(j) = (nRj,y), where 
y being the minimal solution of the equation 



yj{\y-[^jl\)+Fy = E. 



(6.27) 



Moreover, Pp,- = -PR, = -Pq,- 



Equation (6.27) must have two distinct solutions yi^2 {Vi < 1/2) ■ Hi (2/2) corresponds to 



the point from S^{j) with the minimum (maximum) x„. Additionally, since E — Fyi > 
E-Fy2^ Vj {\yi - > Vj (|l/2 - [Ri]„|), we obtain 



■3\n \ ■ 



(6.28) 



To find the maximum of the function ||x — Rj|| on the set S^{j) within the Lagrange 
multipliers method, we introduce the function 



^2{Xn, C, A) = Jc^ + {Xn - + A 



C2 + {Xn - [Rj] J ]+FXn-E 



. (6.29) 



Taking into account inequality (6.28) and the fact that d^2/dc = ^ c = 0, we conclude 



that the maximum of the function ||x — Rj|| on S^{j) is reached at the point pj. 

Let Sj{z) denote a sphere of the radius z centred at Rj, Sj{z) = {x G M" | ||x— Rj|| = z}. 
Consider a sequence of spheres {Sjjrjj + k[Rj — Vj]/^)}Y=o ' where Sj{Rj) = (9supp Vj and 
rjj was introduced in Equation (6.28). Now pick a sequence of points, {'l{k/W)}Y=oy such 
that, ji^k/W) G Sj {rjj + k[Rj — Tij]/W), k = 0, . . . , W. We assume that this sequence is 
a discretization of some differentiable path 7 : [0, 1] — )■ M". According to Equation (6.7), 
the sums, 



w 



Ewh) = v^5^ Jt/,(7(fc/W^))-i5 ||7([A: + l]/W^)-7(^/W^)l|, (6.30) 



k=0 

where we set 7(1 + l/W) = 7(1), obeys the property limvK->oo ^w{l) = -^a(7)- Introduce 
a path: 

g(t) =p,+t[q,-p,]. (6.31) 

Since VA;, g{k/W) G S, {vj + k[R, - r],]/W), [l{k/W)]^ ^ [s{k/W)l and V,{\\g{k/W) - 
R.-ll) = Vjilhik/W) - R,||) U,{j{k/W)) ^ U,{g{k/W))- Moreover, ||7([A; + 1]/W) - 
7(A;/iy)|| ^ \\g{[k + 1]/W) - g{k/W)\\. Therefore, 

Evy(7) ^ Si^(g) ^ LAh) > LA{g). (6.32) 
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Since ^wh) = ^w{g) ^ l{k/W) = g{k/W), k = 0,...,W -1, WW, we conclude that 



path (6.31) is indeed the shortest geodesic that connects the regions S^{j) and 9supp\^-. 
By the same token, the geodesic connecting 9supp Vj and must be a straight hne that 
starts at and ends at Pqj because the potential between these two regions is merely 

To finalize the proof, we shall ''backward propagate" the leading tunnelling trajectory 
starting from the outer turning surface S'^. Let p(x, y) denote the Agmon distance between 
two points for the potential V^(x) = Fxn- Then, it is easy to demonstrate that 



ps(x, Px) = (2/3)V2mF||x - Pxf (g gg) 
The plane T(c) = {x G | a;„ = c} is a surface of constant Agmon distance [Equation 



(6.33)], such that PEiTE/F, S^) = and pe{T{c), S]^) is a strictly increasing function of 



c. Since Hqat — PRatH = Upat — PRatH — Hpiv — qAf|| < \\<ij — PR-jll, Vj 7^ A^, condition 



(6.21) guarantees that increasing c the plane T(c) will "hit" the boundary of suppVA? at 



the point qat. (Note that pe (r(c), S^) = pe {T{c), S^) , E/F < c < [qiv]„.) Moreover, 



the following follows from Equation (6.21) 



{x G M" I [qiv]„ ^ x„, ^ [pn]J H supp Vj = 0, Vj ^ A^, 

which means that the A^*'^ centre is isolated from all the other. Therefore, the shortest 
geodesic must connect the point to the point p^. □ 



Corollary 1. Consider a single electron tunnelling in the potential (6.20), such that as- 
sumptions 1 and 2 of Theorem^ are satisfied, then the leading trajectory is linear (hut may 
not he unique). 

Proof. This corollary follows from the straightforward generalization of the idea of back- 
ward propagation. □ 

Theorem 2. We shall study single electron tunnelling (—00 < E < 0, F > 0) in the 
potential 

K 

^(x) = ^r,(||x-R,||) + Fx.. (6.34) 

Assume that 

1. fj : (0, +00) —7- (—00,0) are differentiahle on (0,+oo) and strictly increasing func- 
tions, such that '^•(0) = —00 and tjl+oo) = 0. 



53 



2. The boundary of the classically forbidden region consists of two disjoin pieces - the 



internal turning surface and the outer one S^. Furthermore, = IJj=i ^eU)' 
^eU) ^ ^eW — ^' 7^ ^> where each S^{j) encircles Rj[^ 

3. B{j)nB{k) = 0, Vj ^ k, andB{])nS^ = 0, Vj, where B{j) = {x e M" | ||x - Rj\\ ^ rj} 
being the ball of radius rj centered at Rj. Here rj = max { ||x — Rj|| | x G S'^(j)} is 
the "radius" of S^Ufl 

Then, the leading tunnelling trajectory (may not be unique) is linear up to a term of 0{X) 
as A — )■ 0, where A = maxj {|7^(Aj)|} and 



A. 



mm 



+1: mill |||R,- - 



2 2 k,k^j 

Here, dj = min | ||x — Rj || | x G 5^} is the Euclidean distance from Rj to S'g. 
Proof. We introduce two auxiliary functions 



(6.35) 







^ X < Aj, 
X ^ A,- 



Vil^{x) 





r,{x) 



^ a; < A 
s ^ Aj. 



J' 



One evidently notices that 



r,{x) = v^^{x)^v^l^{x\ 



(6.36) 



where V^l^ix) is a singular short range potential and yii\x) being a long range tail. The 



purpose of such a partition is to make Vli^ix) satisfy assumption 1 of Theorem 111 and 

'9 



produce Vi^^ix) that obeys the following upper bound: 



Vx. 



We analyze the length of a curve in the Agmon metric [Equation (6.7)]. Since 



E 



K 



|x-R,-||) + Fx„-E + 0(A) 



K 



\ Y.V}l\\\^-R,\\) + Fx^-E + 0{X), 



^ The verb "encircle" should be understood in the following sense. A piece of the inner turning surface, 
^eU) = dCA{j), is a boundary of the classically allowed region, CA{j), associated with centre j, such 
that Rj e CA{j). 

^ The parameter rj can be calculated by means of the method of Lagrange multipliers as it was shown 



in the proof of Theorem [I] [see Equation (6.29)]. 
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under the assumption that A — )■ 0, we have reduced the initial situation to the case of 
single electron tunnelling in the potential 



f/.,(x) = J2 KhHilx - R,||) + Fx„. (6.37) 

i=i 

Let us now utilize assumption 3 to show that 

Aj > Ty (6.38) 

Indeed, on the one hand, B(j) n = ^ ^ dj > r^-; on the other hand, B{j) fl B{k) = 0, 
Vj 7^k, ^ \\Rj - Rfcll - Tk > Tj. 



(6.35) 



assures 



Furthermore, we shall demonstrate that the definition of Aj [Equation (6.35) 
that assumption 2 of Theorem bl for the functions V^l^ (x) holds. According to Equation 



A,- ^ ( II R, - Rfc II - rfc + r,-) /2, j ^ k- 
hence, Aj + A^ ^ ||Rj — Rfc|| ^ suppV^J/J^ flsupp vj^f'' = 0. From Equation (6.35), we also 



obtain Aj ^ dj =^ supp v}p fl 5^ = 0; thus, the outer turning surface for the potential 
Ush{^) should be {x G M" | x„ = E/F}. 

Finally, we have proven the theorem because the potential f/s/i(x) satisfies all the as- 
sumptions of Corollary [Tj □ 

Physical clarifications of Theorems [T] and [2] are due. Assumption 1 of Theorem [T] physi- 
cally implies that Vj are attractive, singular, spherically symmetric short range potentials. 
Assumption 2 of the same theorem requires that the potentials do not merge, i.e., their 
ranges do not overlap. This condition connotes that the classically allowed regions associ- 
ated with the centres Rj [their boundaries are S'^{j)] do not overlap as well. The latter 
statement is proven in Theorem [T] The statement of Corollary [T] can be rephrased as fol- 
lows: leading tunnelling trajectories for a system of non-overlapping, attractive, singular. 



short range potentials are linear. However, if the additional condition (6.21) is satisfied 
then Theorem [T] not only guarantees the uniqueness of the leading tunnelling trajectory, 
but also provides the initial and final points of the trajectory. Assumption 1 of Theorem 
[2] means that fj are attractive, singular, spherically symmetric long range potentials that 
vanish at infinity. Assumptions 2 and 3 of Theorem [2] require the same non-overlapping 
condition for the classically allowed internal regions mentioned above. Physically, Theorem 
[2] says that leading tunnelling trajectories for a system of several such potentials are "al- 
most" linear, and a deviation from being strictly linear is caused by vanishing long tails of 
the potentials; thus, the larger the distance between the centres, the smaller the deviation. 



55 



6.3 The Application of Spherically Symmetric Poten- 
tial Wells to Single-Electron Molecular Tunnelling 



The simplest type of model molecular potentials that allows for full analytical treatment 



is of type (6.20) where 



Cj : < X < Tj, 



X > r 



(6.39) 



J' 



such that Cj < E 



Vj, and it is assumed that S^{j) 



9supp Vj- = {x e M" I ||x|| = rj}. 
(These potentials are not governed by Theorem [l|) Evidently, the leading tunnelling 
trajectories are linear, and moreover, the following equality is valid 



min {pE {qj,Pqj)} 

3 



(6.40) 



where = min^^S'^(j) and pe was defined in Equation (6.33). Let us estimate the 
tunnelling rates within Equation (6.15) for the two dimensional system of two equivalent 



Equations (6.15), (6.33), and (6.40), shows that 



centres of type (6.39) (see Figure 6.1). A straightforward geometrical derivation, using 



rocexp<;-— [Fi?(l 



COS 01 



2E] 



3/2 



(6.41) 



where R is the distance between the potential wells (i.e., the bond length of a model 
molecule) and 6 is the angle between the field and the molecular axis. The obtained 



angular dependent rates are plotted in Figure 6.2 



According to Equation (6.15), rates obtained within the geometrical approach does not 



account for an initial molecular orbital. This technique provides solely the contribution 
from the shape of the barrier, hence, the name - the "geometrical approach." An advantage 
of such a method is that it reduces the calculation of tunnelling rates to a rather simple 
geometrical exercise. 



6.4 Conclusions 

Having introduced the leading tunnelling trajectory as an instanton path that gives the 
highest tunnelling probability, we have proven that leading tunnelling trajectories for multi 
centre short range potentials are linear (Theorem [T| and "almost" linear for multi centre 
long range potentials (Theorem [2]). In a nutshell, these results have been achieved because 
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90° 




270° 



Figure 6.2: The polar plot of the normalized tunnelling rates for the two centre model 
[Equation (6.41 )] as a function of the angle 6. Chosen values of the parameters are F = 0.01 
(a.u.), R = 2 (a.u.), and E = —0.5 (a.u.). 



the multi centre (i.e., molecular) potential is represented as a sum of spherically symmet- 
ric potentials, and such conclusions regarding the shape of the trajectories in the single 
centre (i.e., atomic) case are quite expectable owing to the axial symmetry. An important 
peculiarity of the theorems is that assumptions that they involve are satisfied in major- 
ity situations of current experimental interest. Nevertheless, the proven statements by no 
means exhaust all interesting cases; on the contrary, we have barley scratched the surface, 
and there is plenty of room for further generalizations and expansions. For example, we 
have not discussed the case when internal classically allowed regions associated with the 
centres merge. One may anticipate that the leading tunnelling trajectories still should be 
linear under some additional assumptions (by the argument of the backward propagation 
of the leading tunnelling trajectory). The issue of the uniqueness of the trajectories was 
left nearly undiscussed. 



In any case, one can always employ the fast marching method, discussed in Section [6T 
to obtain numerically some information on the leading tunnelling trajectories. 

The fact that the leading trajectories for long range potentials are not straight lines is 
of vital importance. As in the atomic case 11201 USB [M EHl [ml ^SM USSl Il36] , this 
deviation is crucial for a quantitative treatment [93l [961 EH [HHl 1100] , and sometimes even 
for a qualitative analysis, because it leads to the correct pre-exponential factor of ioniza- 
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tion rates that describes the influence of the Coulomb field of nuclei. However, Theorem 
[2] suggests that the deviations can be accounted for by means of the perturbation theory 
where the zero order approximation being a field free trajectory. This is a part of the 
celebrated PPT approach, widely employed in the literature for analytical calculations of 
atomic Coulomb corrections. Nevertheless, the PPT method requires matching the quasi- 
classical wave function of an electron in the continuum with the bound (field free) atomic 
wave function. This step is a stumbling block for generalization of the PPT approach to 
the molecular case (for the suggestion of a solution to such a problem see Reference |1UU] ). 
Theorem |2] in fact offers a solution to the problem of matching. According to Theorem [2| 
matching should be done on spherical surfaces of radii Aj [Equation (6.35)] centred at the 
nuclei. This is an alternative technique to the method developed in Reference |100j . 

It has been suggested in Reference |170] that molecular photoionization in the tunnelling 
limit may act as a scanning tunnelling microscope (STM). Since rotating a molecule with 
respect to a field direction is analogous to moving the tip of an STM, then the observed 
angular-dependent ionization probability should provide information for a molecule similar 
to the position dependence of the tunnelling current in the STM. We point out that there 
is a resemblance between such a descriptive comparison and our results. As it has been 
shown in Theorem [T] (by the backward propagation of the leading trajectory), the leading 
tunnelling trajectory starts at the atomic centre that is the closest to the barrier exit (i.e., 
the outer turning surface); hence, the qualitative similarity of molecular tunnelling with 
the STM. 
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Chapter 7 
Outlook 



The adiabatic approximation is a convenient, yet intuitive, tool for studying processes 
induced by low frequency laser radiation. In the current work we have studied analytically 
few electron processes such as single-electron ionization and nonsequential double ionization 
(NSDI). 

A future direction, motivated by constantly expanding experimental data, could be 
a general study of many-electron processes. The complexity of the problem decreases 
effectiveness and perhaps even usefulness of any analytical method; ab initio solutions of 
the Schrodinger equation are becoming more and more prominent. Nevertheless, bearing 
a positive attitude towards analytical approaches, we shall discuss two possible directions 
of further applications of the adiabatic approximation. 

First, the adiabatic approximation allows one to reduce the original A^-body problem 
to an effective, or quasi, single body problem within the SFA. We shall demonstrate such 
a reduction on the example of the problem of simultaneous A^ electron ionization [171] . At 
the beginning all the A^ electrons are bounded, i.e.. 



N 



Ef'^ = ~Y.'^^-NI,. (7.1) 



At the end, they are liberated and move under the influence of a linearly polarized laser 
field 

1 ^ 

Ef''\v) = -Y,[^^ + A{^)f. (7.2) 

i=l 



First and foremost. Equation (7.2) is to be simplified. To accomplish that, let us use the 



following simple transformation. For a given set of numbers xi,X2, . . . ,XAr, the standard 
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deviation ax is defined by 



N / N 

i=l \i=l 



xf - Nx^ 



where x = ^ J2f=i denotes the mean; thus, 



TV 



J2^' = N{al + x'). 



(7.3) 



i=l 



Applying Equation (7.3) to Equation (7.2), we obtain 



(7.4) 



Note that the initial (7.1) and final (7.4) energies resemble the initial and final energies 
for single-electron ionization given by Equation (3.1). Let us name the transition from 
Equation (7.2) to Equation (7.4) the mean- deviation parametrization. Having introduced 
the action for simultaneous N electron ionization 



S 



Ne 



^0 



(Ne) 



we can formulate the reduction of such an electron process to the single-electron ioniza- 
tion as 



S'ATe = NS ki^,k±,Ip + 



1 



2 I 2 



(7.5) 



where S is the single-electron action defined in Equation (3.3). Therefore, the ionization 
rates (3.7) are generalized to the case of simultaneous N electron ionization. 

Furthermore, the mean-deviation parametrization can be performed even in the case of 
a "continuous" ensemble of electrons, 

Ef'\v) = l I rf^Pp(p)[p + A(v.)]^ = ^[|^ + A(^)]^ + ^(pi2 + aJ + 4), (7.6) 



where a\\ i are the standard deviations of the momenta and 

II 

given by 

c^^PP(pK±' 



are the mean momenta 



(n = l,2) 



and p(p) is the momentum distribution function of the ensemble. 
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The mean-deviation reduction can be evidently used in many-electron non-sequential 
(correlated) processes and other processes in external fields. In fact, we have already 
applied this trick to the two-electron simultaneous electron emission (SEE) process studied 



in Chapter [5] [see the interpretation of Equation (5.6) . 

However, the role of Coulomb interactions remains the most challenging question in 
many-electron phenomena. The presented mean-deviation parametrization can be most 
naturally introduced only in the SPA, and currently it is not obvious how to account for 
the Coulomb interactions within this formalism. Perhaps, one can develop the SP-EVA-like 
approach to the electron problem; an alternative path could be to use the kinetic theory 

The mean-deviation parametrization, besides being just a computational trick, can be 
used to visualize and process experimental data obtained from many-particle coincidence 
measurements (recently quadruple coincidence measurements have became feasible |173j ). 
According to the mean-deviation parametrization, the averages and standard deviations 
of momenta contain as much information as all the components of the momenta; hence, 
we have a reduction of a many-dimensional data set, which cannot be plotted directly, to 
merely a two-dimensional plot. 

As far as molecular single-electron ionization is concerned, the combination of the 
Hislop-Sigal geometrical approach to tunnelling and the fast marching method (both have 



been presented in Section 6.1) forms a powerful toolkit for studying strong filed phenom- 
ena in polyatomic molecules within the quasiclassical approximation. Purthermore, the 
demonstrated simplicity of the shapes of leading tunnelling trajectories may encourage 
future developments of analytical models of molecular ionization. However, we note that 
the geometrical approach has a fundamental limitation - it does not account for effects 
of molecular orbitals, and there is no a priori way of including these effects. In spite of 
that, one may always attempt to introduce such corrections in a heuristic manner, e.g., 
multiplying the geometrical rates by a Dyson orbital. 

In Chapter [6} we modelled a molecule by a single-electron multi centre potential, hence 
discarding effects of electron-electron correlations. Nevertheless, the geometrical approach 



to tunnelling reviewed in Section |6.1| can account for these effects after an appropriate 
adaptation that is presented in Reference [103J- Intuitively speaking, according to such a 
method, the leading tunnelling trajectory of the system is selected such that the minimum 
number of electrons are displaced during tunnelling. More importantly, the fast marching 
method can be also utilized to obtain this leading tunnelling trajectory. Since correlation 
dynamics of electrons plays an important role in molecular ionization leading to interesting 
novel effects |lU2j . applications of the geometrical ideas developed in Reference [TU3j to 
molecular ionization should be the aim of subsequent studies. 
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Appendix A 

The Phase of Ionization of the 
RecoUiding Electron as a Function of 
the Phase of Recollision 



Employing the SFA, we write the formula corresponding to the digram of NSDI (Figure 



41j) 

\^{t)) ~ fdt, fdtr f d'kU{t,tr)— \kg+) {g+k\ VLih) \gg) x 

expj-^^ [k + A{T)fdT + i\Eg+\{tr-h)+i\Egg\{U-U)^, (A.l) 

where U{t,tr) is the evolution operator of the studied system, is the distance between 
the electrons, Viitb) is the interaction between the ionized electron and the laser field, and 
Eg+ and Egg are energies of the states \g~^) and \gg), respectively. 

We use the saddle point approximation to calculate the integrals over k and tf, in 
Equation (A.l). The phase of the integral over k has the following form: 

Sr{k) = -^-l\k + A{r)fdr. 

The saddle point of this integral is given by 

''^A{T)dr, (A.2) 

h 
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with the restriction tr 7^ tt- Note that generally speaking, k* can be complex since th, as 
will be clarified below, is complex for 7 7^ 0. The phase of the integral over th in Equation 
(A.l) reads 



[k* + A(r)]2rfr + \Eg+\itr - t^) + \Egg\{th - ti). 



Hence, the saddle point tbit^) is a function of tr and given as a solution of the following 
equation 

cos{ujtr) — cos[a;tb(tr)] 



Uj[tr - thitr)] 

where 7 is the Keldysh parameter for the first electron. 



+ sin[u;tfe(tr)] = ±^7, 



7 



E, 



99 J 



2L 



(2) 



(A.3) 



(A.4) 



Methods of computing the saddle points [Equations (A.3) and (A. 7)] have been widely 
discussed in the literature (see, for example. References [63||6ll|65] and references therein). 
We shall use a general and simple approach for identifying correct saddle points between 



different solutions of the saddle-point equations in the complex plane (see Sec. 4.2). 



It is convenient to introduce the following phases: (pi, = ut^ and (^r = ^tr- According 
to Equation (A.3), the saddle point 0f, is a complex double-valued function of (pr which 

^6(+7; ^r) and 0b( 



can be given by 0;,(+7; ^pr) and 0b(— 7; ^r)-, where 7 is the Keldysh parameter (A.4). Here, 
the complex single-valued function 06(7; ^pr 
scendental equation: 



is defined solution of the following tran- 



COS ipr — COS (phi.!] Vr) 

- 4>b{i; Vr) 



+ sin 0b(7; ipr 



(A.5) 



No analytical solution of such an equation is available. The special case of the function 

0b(7; ipr) for 7 = 0, 

— r \- sm (pb{(Pr ) = 0, (A. 6) 

iPr - <Pb[<^r> 

is very important because of the following two reasons. 

First, the function ipbi'Pr) is a real valued function for real ipr (and single-valued for 
any complex ipr)', this allows one to interpret the motion of the first electron in terms of 
classical trajectories. The function (pbifr) is defined on the interval (vr/2, 27r] because only 
during that interval can the free electron recollide with its parent ion. The function ipbifr) 
is bounded in the interval ^ fbiy^r) < 7i"/2. Second, the function ipbifr) can be physically 
understood as a tunneling limit (7 ^ 1) of 0^(7; (fr) (i.e., the low-frequency limit). 
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In terms of the laser phase, the vector potential A{(f) is 

A{{p) = -(F/w) sin 99. 



The difference between Equation (A. 5) and Equation ( A.6[ ), which connect the phase 
of recollision (pr with the phase of ionization (fb{fr) or (j)b{'j','Pr), is also important for the 
last step of NSDI - the release of the two electrons following the recollision at ipr- 



Calculating the integral over tr in Equation (A.l ) by the saddle point approximation, we 
need to obtain the transition point (p^ for negligible 7. It is the solution of the saddle-point 
equation 



such that 



I [k, + A(^°)]^-^ [A(^°) - A((p,(^°))]%jf = 0, (A.7) 



(2) 



where /j 
equation is 



n/2 < Rev?" ^ 27r, 
is the ionization potential of the second electron. For 7 7^ 0, the 



^1) + Hv'r)] ' + ^ N + AivD] ' - \ [A(^°) - A(0,(7; ' + /f = 0, 



where ^6(7; now depends on 7. Note that Equations (A.7) and (A. 8) are basically 
Equation (2.42) written in slightly different notations. 



If the solution of Equation (A.7) on the interval (7r/2, 27r] is real, then direct collisional 
ionization is possible. However, we are interested in the deep quantum regime when the 
following inequality is valid for the second electron: 

> 3.17f/„. 



By introducing the Keldysh parameter for the second electron 72 
write the last inequality as 

72 > 1.26. 



I^^ /{2Up), we can 
(A.9) 



Equation (A.9) physically means that the returning electron does not have enough 
energy to free the second electron. 

When recollision energy is not sufficient for collisional ionization, transition requires 
help from the laser field. Mathematically, the arising integral is similar to those in the 
adiabatic approximation (see Chapter[2]). The energy gap AE{Lpr) in Equations (A.7) and 
(A. 8) plays the role of the transition energy for the non-adiabatic transition [the term 
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Ei(t) — Ef(t) in Equation (2.41)]. The peculiarity of AE{'y,ipr) given by Equation (A. 8) 
is that it need not be real even for real (pr, since in the term [^((^9^) — ^(0fe(7; V'r))] the 
phase (pb is complex. This subtle aspect underscores the important difference between us- 



ing the solutions of Equation (A. 5) or Equation (A. 6) for the phase of ionization 0^. For 
classical trajectories, where i^ti'^r) [Equation (A. 6)] is real for real ipr, AE[ipr) [Equation 



(A. 7)] is also real for real (fr- This is the standard assumption for the adiabatic approxima- 
tion. When the complex phase of ionization (f)i,{'~f](pr) [Equation (A. 5)] is used, i.e., when 



"quantum" trajectories for recoUision are used, AE{'y;(pr) [Equation (A. 8)] need not be 
real for real cpr- In References [TT4j . the Dykhne method (i.e., the adiabatic approximation 
for a system with a discreet spectra) has been generahzed for this case, provided that the 
complex function AE{'j; ip) satisfies the Schwarz reflection principle (recently, this result 
was confirmed and further generalized in Reference [8]). 

The function AE{ip) (A. 7) obeys the Schwarz reflection principle, i.e., AE*{lp*) = 
AE[(p). Hence, we can conclude that if (p^ is the solution that lies in the lower half-plane, 
then {(fr)* is the solution that lies in the upper half-plane. Furthermore, it can be easily 
proven that the following equation takes place for any function AE{(f) which satisfies the 
Schwarz reflection principle and any complex number 

Im / AE{(p)dip = -Im / AE{(f)dip. 

From the previous equation, we can see that the transition points that lie in the lower 
half-plane lead to exponentially large probabilities, which are unphysical. Hereafter, let ip^ 
denote the solution of the equation AE{ip>^) = 0, which is the closest to the real axis and 
lies in the upper-half plane. 

Before continuing our discussion, let us point out the following simple equalities, which 



follow from Equation (A. 5): Re [(j)b{+^] (pr)] = Re [0;,(— 7; (yj^)] and Im [06 (+7; 
— Im [0b(— 7; (pr)] for real ipr- Furthermore, we obtain 



Im / AE{-f;p>)dip 



-Im / AE{ 



-r,^)d^, (A.IO) 



where E{'j; ip) is defined in Equation (A.8). 



Bearing in mind that formula (A.l) must give an exponentially small result (which 
implies that the transition point must be located in the upper-half plane) and taking into 



account Equation (A.IO), we define the phase of ionization in the case of 7 7^ as 



0b(-7; pr) if Im {ip,) > 0, 
Re [0fe(7; p),)] if Im (v^r) = 0, 
(pbl+'j; p>r) if Im {p>r) < 0. 



(A.ll) 



67 



Equation (A. 11) is the most consistent definition of the quantum-mechanical phase of 



ionization of the first electron as a function of the phase of return. Generally speaking, 
there was an ambiguity in selecting the value of $(7;(y9r) for real (pr- However, we have 
chosen it in such a way due to the following reason. The function Im [$(7; ipr)] has a jump 
discontinuity on the real axis, but the function Re [$(7; v?r)] has a removable discontinuity 
that can be removed by employing the equality 

Re [0b(7; v^r)] = ^ [06(7; fr + iO) + 4>b{-T, fr - iO)] (for real (fr)- 

Furthermore, the function $(7; obeys the Schwarz reflection principle [$*(7;(y9*) = 
$(7; fr)], and the following equality takes place: 

$(0; (fr) = fbifr)- 



It is essential that according to definition (A.ll ), the function $(7; fr) is real- valued on the 
real axis and thus allows the identical interpretation in terms of the classical trajectories 



as for (fhlfr)- Therefore, the definition of and inequality (A. 9) are unchanged in the 
case of 7 7^ 0. 
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Appendix B 

Coulomb Corrections in NSDI within 
the Strong-Field Eikonal-Volkov 
Approach 



The key step in dealing with the singularities of the Coulomb potentials during the recol- 
lision is to partition the electron-electron and electron-ion interactions in the two-electron 
Hamiltonian as follows: 



^en — ^en ^en,lng ~\~ ^en,lng 



^ee,lng ^^ee,shr y 
~ Ven,lng ^Ven,shr ■ 



(B.i: 



The potential Vee,ing has a long-range behavior identical to \4e, but no singularity at the 
origin, and ^Vee,shr is singular but short-range potential. The same applies to Ven,ing and 

^Venshr- We choOSC 



AVen,shrir) = Venir)exp{-r/ro), AVee,shriri2) = Veeiri2) exp -ris/r 



(B.2) 



where tq and "^ill be defined later. Note that the partitioning (B.I ) and (B.2) has been 
employed originally in the PPT approach for the problem of single-electron ionization. 

Now, we can write the Hamiltonian as 



Hit) = Hsit) + AVs 



where AK 



shr 



AVee.shr + 



en, shr 



shri 



and is the rest, which includes smoothed Coulomb 



potentials for electron-electron and electron-nuclear interactions. Veering and Ven,ing- 
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To first order in AVghr, the amplitude to find two electrons with momenta ki, k2 at the 
detector at the time t is 



a(ki, k2) = -i 



^ dtr j rf=^k (kiksl Us{t,U)AVshr |^?+k) (k(7+| Usitr^U) \gg) . (B.3) 



Approximations in Equation (B.3) are first order in /SVshr and the assumption that at the 
moment of recolhsion the ion is in its ground state. Both are well justified. 

The next step is to approximate the two parts of the evolution: before tr and after tr- 
The key component for correlated spectra is the second part - after tr. The main aspect 
of the first part of the evolution - prior to t^. - is to supply an active electron with the 
required energy. 

To simplify the amplitude (kik2| tJs{t,tr)AVshr \g^^)i we insert the decomposition of 
unity, 

6(ki,k2,k,tr) = 



(kik2|f/.(t,t.)AV;,,, |(7+k) 
j j dh^dh^ (kiksl Us{t,t,) Irirs) (rir2| AK/^r j^^+k) 

j j dhidh2 (ki + A(t,), ks + A(t,)| rir2) (rir2| AVshr |^+k) 



X exp 



—I 



^[ki + A(r)]2 + i[k2 + A(r)]2+ 



+V;e,/ns(ri2(r)) + Ven,lng{ri{T)) + Kn,«ng (r2 (r) ) | rfrj . (B.4) 

Here we have applied the strong-field eikonal-Volkov approach |133j (SF-EVA) method 
|133j . The integral from the nonsingular parts of the electron-electron and electron-ion 
interactions are calculated along the trajectories in the laser field only. The trajectories 



ri,2(t) = ri,2+ / [ki,2 + A(r)]rfr 



(B.5) 



and ri2(t) = ri(t) — r2(^) begin at the positions ri,r2 at instant tr- The bra-vectors 
(ki_2 + A(tr)| are plane waves. Their distortion by the electron-electron and electron-core 
interactions appears in the (ri, r2)-dependent exponential phase factors in Equation (B.4). 

Since AVshr is a short-range potential and \g^) is limited within a characteristic ionic 
radius, the term (rir2| AVghr \9^'^) allows us to fix the initial values of ri and Y2. The 
characteristic radius for the partitioning of the Coulomb potential into the short-range and 



long-range parts is set as tq 



,(0) 

12 



1/ I £'3+ 1, where Eg+ denotes the energy level of the 
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second (b ound ) electron. Therefore, we pull the exponential factor out of the integral in 



Equation (B.4) with ro = r 



(0) 
12 



1/ \Eg+ \ and ri = r2 = 0, 



-zy^'|i[ki + A(r)]^+ 

(B.6) 



6(ki,k2,k,t,) ^ (kl + A(^,)k2 + A(^,)|AK/^r|^+k)exp 



+ -[k2 + A(r)] + V;e,«ng(ri2(i")) + V;„,i„g(ri(r)) + V;„,i„(,(r2(r)) | rfr 

Effects of the long-range tails of V^e and Ven appear in the exponent while the collisional 
transition is govered by the short-range interaction |133] . The states |k), so far, represent 
any convenient basis set of continuum sates in the laser field. 

To simphfy the amplitude 

c(k,tr) = {kg+\ Us{tr,ti) \gg) , 

we note that the second electron is bound during the whole evolution, and hence we can 
simplify c(k, tj) using single active electron approximation. In this approximation, Us{tr, ti) 
describes one-electron dynamics in the self-consistent potential of the ionic core, 

Kc(ri) = {g^ \ Vee,«ng(ri2) + Ven,lng{ri) + Ven,lngi^2) \g^) ■ 

The effective Hamiltonian for evolution between tj and tr is 

Hscirut) = Ki + Kc(ri) + ^^(ri,^), 

where Ki is the kinetic energy operator and Vi(ri, t) is the interaction with the laser field. 
Now the amplitude c(k, tj.) becomes 



tr 



dtb (k| Uscitr, tb)VL{ri,tb) Igo) exp [i \Egg\ {U - U)] 



(B.7) 



where {go) = {gt loo) is proportional to the Dyson orbital between the ground states of 
the neutral and ion and Egg is the energy of the ground state l^'^'). 



The ionic potential contributes to the propagator in Equation (B.7) twice: when the 
electron leaves the atom near tf, and when it returns to the ionic core near tr- The contribu- 
tion "on the way out" introduces standard Coulomb correction |1191I12U[I1211I122"| 11351 [TM] 
to the ionization amplitude and hence affects the overall height of the final two-electron 
distribution. The contribution of Vgc "on the way in" affects the spatial structure of the 
recoUiding wave packet. As shown in References |175| 1133] , for short collision times the 
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Coulomb-laser coupling is small and Vsc "on the way in" can be included in the adiabatic 
approximation, 



c{k,tr)\K^g^) ^ -iRc / dtb\ke^g^) {k + A{U) - A{tr)\VL{tb)\9D) 



exp 



[k + A(r) - A{tr)]^dr + i\Eg+ \ {U - h) + i\Egg\ {U - U 



(B.8) 



Here |ket,) is the field-free continuum wave function in the eikonal approximation, which in- 
cludes distortions of the incoming plane wave with asymptotic momentum k, (k -|- A{tb) — A(i 
is a plane wave, and Rc is the Coulomb correction to the ionization amplitude which com- 
pensates for approximating (k + A{tb) — A(tj.)| with a plane wave in the matrix element 
{k + A{tb) - A{U)\VLih)\gD). 



Now, putting together Equations (B.8) and (B.6) and changing the integration variable 
k — )■ k + A{tr), we arrive at 



a(ki,k2 



dtr / dtb / d^k / d vid r2 



X exp 



-- / [k + A{T)]^dT + i\Eg+\{tr-tb)+l\Egg\{tb-ti 

i[ki + A(r)]2 + i[k2 + A(r)]2}dr 

-i {V;e,/ns(ri2(l")) + V;n,«n(,(ri(r)) + Kn,«ng (r2 (^) ) } dr 

X (ki + A(t,), k2 + A(t,)|rir2) 

X (nral /W,hr |^+, + A{U)) Rc (k + A{tb)\ Viih) |^d) ■ 



(B.9) 



Note that if the Coulomb corrections Vee,ing and Kn,«ng are ignored in the exponent of 
Equation (B.9), then Equation (B.9) coincides with Equation (A.l) within exponential 
accuracy. 
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Appendix C 



Upper Bounds for Matrix Elements 
and Transition Amplitudes 



We derive a multi-dimensional generalization of the Landau method of calculating qua- 



siclassical matrix elements [Equation (C.4)], and we also estimate perturbation theory 



transition amplitudes [Equations (C.6) and (C.9)] in terms of the Agmon distance. 

For the sake of simplicity, the argument x will be omitted in some equations below. 
Throughout this Appendix, we assume that the Agmon upper bounds [150J for bound 
states {ipn) are valid, i.e., Ve > 3cn = c„(e), < c„ < cxd, such that 



id] 



where p„ = pe„. 

Let us choose an arbitrary e > 0. Employing the Schwartz inequality and assumption 



(C.l ), we obtain 



{l-e){pp + pg)^*^^^-{l-e){pp+Pq)Y^^ 



where e > e' > and 



B: 



|^|'g-'(l-e)(Pp + P<7)^-j 



(C.2) 



(C.3) 
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The integral / exp[— 2(e — e')(pp + Pq)]c?x converges for all p and q. Therefore, we have 
proven that Ve > 3c = c(e), < c < oo, such that 



^ cB. 



p,q^ 



(C.4) 



which is the same as Equation (6.12). 



Now let us study the problem of estimating of transition amplitudes defined by means 
of the time dependent perturbation theory. Hereinafter, we assume that a quantum system 
under scrutiny has no continuum spectrum, and we shall manipulate with all the series and 
integrals heuristically - assuming that they all converge, or alternatively, assuming that 
they are over finite range. We illustrate our idea by estimating the second order amplitude 
since generalization to higher orders will be evident. 

The second order transition amplitude within the time dependent perturbation theory 
reads 



dt I dt' I rfx(ixVfi„(x')e~*^-'"(*^-*') 



-iEi„{t-ti) 



xV{x.')K{x.'t'\xt)V{x.)^lj^n{x.)e 



(C.5) 



where all the ip's are eigenstates of the system and K is the propagator, which can be 
written as 



ir(xY|xt) = 5^^„(x')C(x)e- 



-iEnit'-t). 



whence, 



|ir(xY|xt)| ^5^|V^„(x')^„(x)|. 



Using such a simple estimate as well as inequality (C.l), we obtain 



(C.6) 



where M = CmCjin max„ {c^} /2, < M < oo. 

However, there is no need to confine ourself to the case when the initial and final states 
are eigenstates. The same idea applies to the general case of the initial {(pin) and final 
{(f) fin) states being represented as linear expansions in the basis of the bound eigenstates. 



(C.7) 
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Let us found an upper bound for the first order transition amplitude, which is as follows 

n,n' 

xl^(x)(V^„, ^l^^,{-K.)e-^''r.'it-u\ (C.8) 

Whence, we readily obtain 

N(i)| 

j—^ ^ Y,CnCn'\{(t>finm{i^n'\<t>in)\Bln'^MY,Bl^,, (C.9) 
n,n n,n 

where M = max„^„/ {c„Cn/|(0/i„ |V^„) (V',,/ |}, < M < oo. 
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